Lorentz and CPT violation in neutrinos 
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A general formalism is presented for violations of Lorentz and CPT symmetry in the neutrino 
sector. The effective hamiltonian for neutrino propagation in the presence of Lorentz and CPT 
violation is derived, and its properties are studied. Possible definitive signals in existing and future 
neutrino-oscillation experiments are discussed. Among the predictions are direction-dependent ef- 
fects, including neutrino-antineutrino mixing, sidereal and annual variations, and compass asymme- 
tries. Other consequences of Lorentz and CPT violation involve unconventional energy dependences 
in oscillation lengths and mixing angles. A variety of simple models both with and without neutrino 
masses are developed to illustrate key physical effects. The attainable sensitivities to coefficients for 
Lorentz violation in the Standard-Model Extension are estimated for various types of experiments. 
Many experiments have potential sensitivity to Planck-suppressed effects, comparable to the best 
tests in other sectors. The lack of existing experimental constraints, the wide range of available 
coefficient space, and the variety of novel effects imply that some or perhaps even all of the existing 
data on neutrino oscillations might be due to Lorentz and CPT violation. 
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I. INTRODUCTION 

The minimal Standard Model (SM) of particle physics 
offers a successful description of most processes in Na- 
ture but leaves unresolved several experimental and the- 
oretical issues. On the experimental front, observations 
of neutrino oscillations have accumulated convincing ev- 
idence that the description of physical properties of neu- 
trinos requires modification of the neutrino sector in the 
minimal SM. Most experimental results to date can be 
described theoretically by adding neutrino masses to the 
minimal SM, but a complete understanding of the ex- 
isting data awaits further experimentation. On the the- 
oretical front, the SM is expected to be the low-energy 
limit of a more fundamental theory that unifies quan- 
tum physics and gravity at the Planck scale, mp ~ 10^^ 
GeV. Direct measurements at this energy scale are im- 
practical, but suppressed low-energy signatures from the 
anticipated new physics might be detectable in sensitive 
existing experiments. 

In this work, we address both these topics by study- 
ing effects on the neutrino sector of relativity violations, 
a promising class of Planck-scale signals. These viola- 
tions might arise through the breaking of Lorentz sym- 
metry and perhaps also the breaking of CPT symmetry 
[1]. Since the SM is known to provide a successful de- 
scription of most physics at low energies compared to 
the Planck scale, any such signals must appear at low 
energies in the form of an effective quantum field the- 
ory containing the SM. The general effective quantum 
field theory constructed from the SM and allowing arbi- 
trary coordinate-independent Lorentz violation is called 
the Standard-Model Extension (SME) [2]. It provides a 
link to the Planck scale through operators of nonrenor- 
malizable dimension [3, 4]. Since CPT violation implies 
Lorentz violation [5], this theory also allows for general 
CPT breaking. The SME therefore provides a realistic 



theoretical basis for studies of Lorentz violation, with or 
without CPT breaking. 

The lagrangian of the SME consists of the usual SM 
lagrangian supplemented by all possible terms that can 
be constructed with SM fields and that introduce viola- 
tions of Lorentz symmetry. The additional terms have 
the form of Lorentz-violating operators coupled to coef- 
ficients with Lorentz indices, and they could arise in a 
variety of ways. One generic and elegant mechanism is 
spontaneous Lorentz violation, proposed first in string 
theory and field theories with gravity [6] and then gen- 
eralized to include CPT violation [7]. Another popular 
framework for Lorentz violation is noncommutative field 
theory, in which realistic models form a subset of the 
SME involving operators of nonrenormalizable dimension 
[8]. Other proposed sources of Lorentz and CPT viola- 
tion include various non-string approaches to quantum 
gravity [9], random dynamics [10] and multiverses [11]. 
Planck-scale sensitivity to the coefficients for Lorentz vi- 
olation in the SME has been achieved in various exper- 
iments, including ones with mesons [3, 12, 13], baryons 
[14-16], electrons [17, 18], photons [19-22], and muons 
[23]. However, no experiments to date have measured 
neutrino-sector coefficients for Lorentz violation. 

Here, we explore neutrino behavior in the presence of 
Lorentz and CPT violation using the SME framework. 
The original proposal for Lorentz and CPT violation in 
neutrinos [2] has since been followed by several theoret- 
ical investigations within the context of the SME [24- 
29], most of which have chosen to restrict attention to 
a small number of coefficients. A comprehensive the- 
oretical study of Lorentz and CPT violation in neutri- 
nos has been lacking. The present work partially fills 
this gap by applying the ideas of the SME to a gen- 
eral neutrino sector with all possible couplings of left- 
and right-handed neutrinos and with sterile neutrinos. 
We concentrate mostly on Lorentz-violating operators 
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of renormalizablc dimension, which dominate the low- 
energy physics in typical theories, but some generic con- 
sequences of Lorcntz-violating operators of nonrenormal- 
izable dimension are also considered [3, 4, 30] . The effec- 
tive hamiltonian describing free neutrino propagation is 
obtained, and its implications are studied. The formalism 
presented in this work thereby provides a general theo- 
retical basis for future studies of Lorentz and CPT viola- 
tion in neutrinos. We also illustrate various key physical 
ideas of Lorentz and CPT violation through simple mod- 
els, and we discuss experimental signals. Our primary 
focus here is on oscillation data [31], but the formalism 
is applicable also to other types of experiments including 
direct mass searches [32] , neutrinoless double-beta decay 
[33] , and supernova neutrinos [34] . 

Several features of Lorentz and CPT violation that we 
uncover are common to other sectors of the SME, in- 
cluding unconventional energy dependence and depen- 
dence on the direction of propagation. We also find 
that Lorentz- violating neutrino-antineutrino mixing with 
lepton-number violation naturally arises from Majorana- 
like couplings. These features lead to several unique sig- 
nals for Lorentz and CPT violation. For example, the 
direction dependence potentially generates sidereal vari- 
ations in terrestrial experiments as the Earth rotates, an- 
nual variations in solar-neutrino properties, and intrinsic 
differences in neutrino flux from different points on the 
compass or different angular heights at the location of the 
detector. The unconventional energy dependence pro- 
duces a variety of interesting potential signals, including 
resonances in the vacuum [25, 29] as well as the usual 
MSW resonances in matter [35]. 

Experiments producing evidence for neutrino oscilla- 
tions to date include atmospheric-neutrino experiments 
[36], solar-neutrino experiments [37-42], reactor experi- 
ments [43], and accelerator-based experiments [44, 45]. 
Most ciirrent data are consistent with the introduction 
of three massive- neutrino states, usually attributed to 
GUT-scale physics. However, as we demonstrate in this 
work, the possibility remains that the observed neutrino 
oscillations may be due at least in part and conceivably 
even entirely to Lorentz and CPT violation from the 
Planck scale. In any event, experiments designed to test 
neutrino mass are also well suited for tests of Lorentz and 
CPT invariance, and they have the potential to produce 
the first measurements of violations of these fundamental 
symmetries, signaling possible Planck-scale physics. 

The organization of this paper is as follows. Section 
II presents the basic theory and definitions, obtaining 
the effective hamiltonian for neutrino propagation and 
discussing its properties. Issues of experimental sensitiv- 
ities and possible constraints from experiments in other 
sectors are considered in Section III. Certain key features 
of neutrino behavior in the presence of Lorentz and CPT 
violation are illustrated in the sample models of Section 
IV. Some remarks about both generic and experiment- 
specific predictions are provided in Section V. Through- 
out, we follow the notation and conventions of Refs. [2, 4]. 



II. THEORY 
A. Basics 

Our starting point is a general theory describing N 
neutrino species. The theory is assumed to include all 
possible Majorana- and Dirac-type couplings of left- and 
right-handed neutrinos, including Lorentz- and CPT- 
violating ones. The neutrino sector of the minimal SME 
is therefore included, along with other terms such as 
those involving right-handed neutrinos. 

We denote the neutrino fields by the set of Dirac 
spinors {ve, ^t,- ■ ■} and their charge conjugates by 
{v^c = , v^c = , Vt-c = ,. . . }, where charge con- 
jugation of a Dirac spinor is defined as usual: = Cd"^ . 
By definition, active neutrinos are detected via weak in- 
teractions with left-handed components of {i/g, i^fj., Vt}- 
Complications may arise in the full SME, where Lorentz- 
violating terms alter these interactions and can modify 
the detection process. However, such modifications are 
expected to be tiny and well beyond the sensitivity of cur- 
rent experiments. In contrast, propagation effects can be- 
come appreciable for large baselines. We therefore focus 
in this work on solutions to the Lorentz-violating equa- 
tions of motion that describe free propagation of the 
neutrino species. 

It is convenient to place all the fields and their conju- 
gates into a single object va, where the index A ranges 
over the 2N possibilities {e, /u, r, . . . , e*-^, /z'", t*-^, . . .}. 
This setup allows us to write the equations of motion 
in a form analogous to the Lorentz-violating QED exten- 
sion [2, 4], and it can readily accommodate Dirac, Ma- 
jorana, or more general types of neutrinos. Our explicit 
analysis in this section is performed under the assump- 
tion that Lorentz-violating operators of renormalizablc 
dimension dominate the low-energy physics. Then, the 
general equations of motion for free propagation can be 
written as a first-order differential operator acting on the 
object VA'- 

{ir^ABd. - Mab)vb = 0. (1) 

Here, each constant quantity T^^, Mab is also a 4 x 4 
matrix in spinor space. Note that the usual equations 
of motion for Dirac and Majorana neutrinos are special 
cases of this equation. 

The matrices and Mab can be decomposed using 
the basis of 7 matrices. We define 

Tab = l^^AB+cf^ABlti + d^ABl^ln 

+eAB + ifABlh + ^ABf^Xl^^ 

Mab = mAB + irn^ABlb 

+"-ABlf^ + bAB^5l,, + ^H'^b'^^^. (2) 

In these equations, the masses m and ms are Lorentz 
and CPT conserving. The coefficients c, d, H are CPT 

conserving but Lorentz violating, while a, 6, e, /, g are 
both CPT and Lorentz violating. Requiring hermiticity 
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of the theory imposes the conditions FJ^^ = j'^i^Bj^YY' 
and Mab = 7°(-Mba)^7°) which impHes aU coefficients 
arc hcrmitian in generation space. 

The above construction carries some redundancies that 
stem from the interdependence of v and f'-^ . This imphes 
certain symmetries for and M. Note first that charge 
conjugation can be written as a Unear transformation 
on va'- i^A— ^ABi^B, where C is the symmetric matrix 
with nonzero elements C^^c = C^^c = C^^c = • • • = 
1. Then, in terms of C and the spinor matrix C, the 
interdependence of u and impHes the relations 

Tab = -CacCbdC{T'^q)'^C~^ , 
Mae = CAcCBDC{MncfC-\ (3) 

where the transpose T acts in spinor space. Suppressing 
generation indices, this translates to 

cA"' = C(c^'^)^C, m = C(m)^C, 

= -Cie^fC, a'' = -Cia^fC, (4) 
= -cId^C, = C(6^)^C, 

where now the transpose T acts in generation space. 
Note that the overall sign in the above equations are 
chosen to match their derivation within the conventional 
lagrangian formalism involving anticommuting fermion 
fields. 

Equation (1) provides a basis for a general Lorentz- 
and CPT-violating relativistic quantum mechanics of 
freely propagating neutrinos. However, the unconven- 
tional time-derivative term complicates the construction 
of the corresponding hamiltonian. This difficulty also 
arises in the minimal QED extension, but it may be over- 
come [4] if there exists a nonsingular matrix A satisfying 
the relationship A^'y'^F'^A = 1. The field redefinition 
va = AabXb then allows the equations of motion (1) to 
be written as {iSABdo — 'Hab)xb = 0, where the hamil- 
tonian is given by H = -A^^°{ir^dj - M)A. 

Denoting SV^ and SM as the Lorentz-violating por- 
tions of and M, and under the reasonable assumption 
that |(5r"| < 1, a satisfactory field redefinition is given by 
the power scries ^ = (l+j'^ST^)-^/^ = 1 - ^-f^ST^ + ■ ■ ■ . 
Separating the hamiltonian H into a Lorentz-conserving 
part Ho and a Lorentz-violating part 6H, which we as- 
sume is small relative to TCq, wc can use the above ex- 
pression for A to obtain an expansion of STi in terms 
of Ho and coefficients for Lorentz violation. Explicitly, 
at leading order in coefficients for Lorentz violation, we 
obtain 

SH = -^{j%T'>Ho+Hol°6T°)--f°{i6T^dj - 6M). (5) 

This expression is therefore the basis for a general study 
of leading-order Lorentz and CPT violation in the neu- 
trino sector. 

At this stage, prior to beginning our study of Eq. 
(5), it is useful to review the properties of the Lorentz- 
conserving hamiltonian [46, 47] 

Ho = -i\i^'dj - Mo). (6) 



The Lorentz-conserving dynamics is completely deter- 
mined by the mass matrix Mg, which in its general form 
can be written 

Mo = m + zm575 = ttilPl + mnPR , (7) 

with ruR = {niLy ^ m + im^ and Pl = \{l - 75), Pr = 
^(1 -|- 75). The components of the matrix tur = m\ 
can be identified with Dirac- or Majorana-type masses 
by separating tur into four N x N submatrices. It is 
often encountered in the form of the symmetric matrix 

The matrices R and L are the right- and left-handed 

Majorana-mass matrices, while D is the Dirac-mass ma- 
trix. In general, R, L and D are complex matrices re- 
stricted only by the requirement that R and L are sym- 
metric. Note that a left-handed Majorana coupling is 
incompatible with electroweak-gauge invariance. In con- 
trast, Dirac and right-handed Majorana couplings can 
preserve the usual gauge invariance. 

It is always possible to find a basis in which the mass 
matrix Mq is diagonal. Labeling the fields in this basis 
by XA', where A' = 1, . . . , 2N, then the unitary transfor- 
mation relating the two bases can be written as 

Ua'A = Va'aPl + {VCYa,aPr, (9) 

where is a 2N x 27V unitary matrix. Here, it is un- 
derstood that Ua'A carries spinor indices that have been 
suppressed. In the new basis, the mass matrix tola'B' = 
rriRA'B' — 'n^(A')^A'B' is diagonal with real nonnega- 
tive entries. The neutrinos xa' = Xa' — Va'aPlXa + 
Va'A^RXa Majorana particles, regardless of the form 
of Mq. 

B. Effective hamiltonian 

The discussion above applies to an arbitrary number of 
neutrino species and an arbitrary mass spectrum. Since 
a general treatment is rather cumbersome, we restrict at- 
tention in what follows to the minimal physically reason- 
able extension with N = 3. For definiteness, we also as- 
sume a standard seesaw mechanism [48] with the compo- 
nents of R much larger than those of D or L. This mech- 
anism suppresses the propagation of right-handed neu- 
trinos, so the analysis below also contains other Lorentz- 
and CPT-violating scenarios dominated by light or mass- 
less left-handed neutrinos, including the minimal SME. 

Ordering the masses m(^A') from smallest to largest, 
we assume that m(i), m(2), m(3) are small compared to 
the neutrino energies and possibly zero, and that the re- 
maining masses m(4) , m(5) , m(6) are large with the corre- 
sponding energy eigenstates kinematically forbidden. In 
this situation the submatrix Va'a, where a = e, fj,,T and 
a' = 1,2,3, is approximately unitary. 
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To aid in solving the equations of motion, we define 

XAit;x) = J ^XA{t;p)e'P^-^ , 

XA{t;p) = bA{t;p)uL{p) + {Cd)A{t;p)uii{p) 

+(C6)l(t; -p)vn{-p) + d\{t; -p)vl{-p). 

(10) 

This is chosen to satisfy explicitly the charge- 
conjugation condition — CabXb- The spinor basis 
{uhiv), Uii{p), Vfi{~p),VL{—p)} obeys the usual relations 
for massless fermions, with vr^l{p) = CvJ^ j^{p). It has 
eigenvalues of the helicity operator 757*^7 • p/\'p\ given 
by {—,+,—,+} and eigenvalues of the chirality operator 
75 given by {—,+,+,—}. For simplicity, we normalize 
with ujjM/g = vl^vp = Sa/s for a,0 = L. R. The definition 
(10) implies that the amplitudes be,^,T may be approxi- 
mately identified with active neutrinos and de,jx,T with ac- 
tive antineutrinos. The remaining amplitudes b^c^^c ^.c 
and deC ^^C T-c cover the space of sterile right-handed neu- 
trinos, but a simple identification with flavor neutrinos 
and antineutrinos would be inappropriate in view of their 
large mass. 

In the mass-diagonal Majorana basis, wc restrict at- 
tention to the propagating states consisting of the light 
neutrinos. Taking the hamiltonian in this basis, 

'Ha'b'ip) = l^in ■P + m(a'))5a'b' +5'Ha'b'{p), (H) 

and applying it to Xb'{t;p) = Ub'BXB{t;p} yields the 



equations of motion in terms of the amplitudes b and 
d. The result takes the form of the matrix equation 

[iSa^b'do-H^'b'ip)] =0, (12) 

\di{t;-p)J 

where for convenience we have defined bb' = Vb'BbB and 
db' = VygdB, and where Ha'b' is the spinor-decomposed 
form of Ha'b'. 

The propagation of kinematically allowed states is 
completely determined by the amplitudes and da'- 
However, for purposes of comparison with experiment 
it is convenient to express the result using the ampli- 
tudes associated with active neutrinos, 6e,/^,r and de.f^.j- 
The relevant calculation is somewhat lengthy and is de- 
ferred to Appendix A. It assumes that the submatrix 
Va'a is unitary, and it neglects terms that enter as small 
masses m^a') multiplied by coefficients for Lorentz viola- 
tion, since these are typically suppressed. The calcula- 
tion reveals that the time evolution of the active- neutrino 
amplitudes is given by the equation 

(*?.;g)--HM.).(s:<S:g), (13, 

where /loir is the effective hamiltonian describing flavor 
neutrino propagation. To leading order, it is given by 



{Ks)ab = Wab 1+^ 



1 \ 1 A ifh%b 



2[p\ \ (m 
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Jab 



[iaLrPt.-icLyPf.P.]ab -iV2p^,{e+U{g^''"'p,-HnCU\ 
1^ I /o ^ ^*^/'«u^'c^„ I zJ"u^'^/^^* ^/'«_^w^ f^_\Liv^ ^ I' v / 



\p\ V iV2p^{e+)l[{g^'^^p. + H^CMb - {cLrv^^vAlb 



where we have deflned {cL)ab = (c + d)'^'^ and (aL)^^ = 
(a + 6)^^ for reasons explained below. The approximate 
four momentum may be taken as = (|p|; —p) at 
leading order. The Lorentz-conserving mass term re- 
sults from the usual seesaw mechanism with = m;mj', 
where m; is the light-mass matrix mi = L — DR~^D . 
The complex vector (e+)^ satisfles the conditions 

p^{e+r -p''{e+Y = ie^-'"'Ppie+U 

{e+ne+): = -1. (15) 

A suitable choice is {e+Y = -^(0; ei + ^62), where ei, €2 
are real and ei, 62} form a right-handed orthonor- 

mal triad. Note that {e+Y and {€-)" = {e+Y* is analo- 
gous to the usual photon helicity basis. The appearance 
of these vectors reflects the near-definite helicity of active 
neutrinos. The vectors ei and €2 can be arbitrarily set 



by rotations or equivalently by multiplying {e-^-)'^' by a 
phase, which turns out to be equivalent to changing the 
relative phase between the basis spinors ul and ur. 

Only the diagonal kinetic term in hcs arises in the 
minimal SM. The term involving {fn^)ab encompasses 
the usual massive-neutrino case without sterile neutri- 
nos. The leading-order Lorentz- violating contributions to 
neutrino-neutrino mixing are controlled by the coefficient 
combinations (a + &)af, and (c+d)'^'^. These combinations 
conserve the usual SU(3)xSU(2)xU(l) gauge symmetry 
and correspond to the coefficients (aL)^^ and (cl)^^ in 
the minimal SME. Note that the orthogonal combina- 
tions (a — 6)^;, and (c — d^^ also conserve the usual 
gauge symmetry, but they correspond to self- couplings 
of right-handed neutrinos and are therefore irrelevant 
for leading-order processes involving active neutrinos. 
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The remaining coefScicnts, {g^'^'^C)ab and (H'^'^C)ab, ap- 
pear in /icir through Majorana-hke couphngs that violate 
SU(3) X SU(2) X U(l) gauge invariance and lepton-number 
conservation. They generate Lorentz-violating neutrino- 
antineutrino mixing. 

Some combinations of coefficients may be unobserv- 
able, either due to symmetries or because they can be 
removed through field redefinitions [2, 4, 49, 50]. For 
example, the trace component rj^vicL)^^ is Lorcntz in- 
variant and can be absorbed into the usual kinetic term, 
so it may be assumed zero for convenience. In fact, even 
if this combination is initially nonzero, it remains ab- 
sent from the leading-order effective hamiltonian because 
the trace of p^Pv vanishes. Other examples of unob- 
servable coefficients include certain combinations of g^'^'^ 
and H^^" . The antisymmetry properties g^'^'^ = —g'^^'^, 
H^^'^ = — iJ'^^ and the properties of (e+)^ can be com- 
bined to prove that the physically significant combina- 
tions of g^"^"^ and H^'^ are given by the relations 

P^{e+).H'^-' = b1(e+).i^^ (16) 
where we have defined 

jjt^ = ijO"^ + i eO^^piJw . (17) 

Only these combinations appear in hes and are relevant 
to neutrino oscillations. 

In deriving Eq. (14), we have focused on operators of 
renormalizable dimension, which involve linear deriva- 
tives in the equations of motion and a single power of 
momentum in the hamiltonian. Operators of nonrenor- 
malizable mass dimension n > 4 are also of potential im- 
portance [3, 4]. They appear as higher-derivative terms 
in the action, along with corresponding complications in 
the equations of motion and in the construction of the 
hamiltonian. An operator of dimension n is associated 
with a term in the action involving d — n — S derivatives, 
and the associated terms in the effective hamiltonian in- 
volve d powers of the momentum. The corresponding 
coefficient for Lorentz violation carries d + 2 or fewer 
Lorentz indices, depending on the spinor structure of the 
coupling and the number of momentum contractions oc- 
curring. For the case n > 4, we gencrically denote the 
coefficients by (kd)^'"- These coefficients have mass di- 
mension 1 — d. Note that, depending on the theory con- 
sidered, the mechanism for Lorentz and CPT violation 
can cause them to be suppressed by d-dependent powers 
of the Planck scale [3, 4] . Some effects of operators with 
d = 2 have been considered in the context of quantum 
gravity in Ref. [30] . 

The mixing described by Eq. (14) or its generalization 
to operators of dimension n > 4 can be strongly energy 
dependent. For example, any nonzero mass-squared dif- 
ferences dominate the hamiltonian at some low-energy 
scale. However, while mass effects decrease with energy. 



Lorentz-violating effects involving operators of renormal- 
izable dimension remain constant or grow linearly with 
energy E and so always dominate at high energies. For 
instance, the contributions from a mass of 0.1 eV and a 
dimensionless coefficient of 10"^'' are roughly compara- 
ble at an energy determined by E'^ ~ (0.1 eV)^/(10~"'^^), 
or ~ 30 MeV. Below this energy the mass term dom- 
inates, while above it the Lorentz-violating term does. 
Similarly, a dimension-one coefficient of lO^"'^'"' GeV has 
a transition energy ~ 10 keV. More generally, effects 
controlled by the coefficients (fcd)^ '' for Lorentz violation 
involving operators of dimension n = d + 3 grow as E"^. 

Although the perturbative diagonalization leading to 
Eq. (14) is valid for dimensionless coefficients much 
smaller than one and for energies much greater than any 
masses or coefficients of dimension one, at sufficiently 
high energies issues of stability and causality may require 
the inclusion of Lorentz-violating terms of nonrenormal- 
izable dimension in the theory. In the context of the 
single-fermion QED extension, for example, a dimension- 
less c°° coefficient can lead to issues with causality and 
stability at energies ~ mfermion/\/cO^ unless the effects 
of operators of nonrcnormalizable dimension are incor- 
porated [4]. A complete resolution of this issue would 
be of interest but lies beyond our present scope. It is 
likely to depend on the underlying mechanisms leading 
to mass and Lorentz violation, and it may be complicated 
further by the presence of multiple generations and the 
sterile neutrino sector. We limit our remarks here to 
noting that the values of the coefficients for Lorentz vio- 
lation considered in all the models in this work are suffi- 
ciently small that issues of stability and causality can be 
arranged to appear only beyond experimentally relevant 
energies. In any case, the renormalizable sector provides 
a solid foundation for the basic treatment of Lorentz and 
CPT violation in neutrinos. 

C. Neutrinos in matter 

In many situations, neutrinos traverse a significant 
volume of ordinary matter before detection. The re- 
sulting forward scattering with electrons, protons, and 
neutrons can have dramatic consequences on neutrino 
oscillations [51]. These matter interactions can read- 
ily be incorporated into our general formalism. Since 
the effective lagrangian in normal matter is given by 

A£mattor = -V2G prieVel^ PlI^c + [G / V^)!^,^'^ Pll^a, 

matter effects are equivalent to contributions from CPT- 
odd coefficients 

(aL,eff)ee = GF(2ne - n„)/\/2, 

(aL,eff)°^ = {aL,ef£)rr = -GFnn/V2, (18) 

where Ue and n„ are the number densities of electrons 
and neutrons. Adding these terms to the effective hamil- 
tonian (14) therefore incorporates the effects of matter. 

For some of the analyses of Lorentz violation below, it 
is useful to review the treatment of matter effects in solar 
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and atmospheric neutrinos. Consider first solar neutri- 
nos. These are produced in several processes that gen- 
erate distinct, well-understood Ve spectra. The most no- 
table are the pp spectrum with a maximum energy of 
about 0.4 MeV, and the ®B spectrum with a maximum 
of about 16 McV [52]. For Va I'b mixing scenarios, 
the contribution from n„ is the same for all species and 
therefore can be ignored. However, n„ may be impor- 
tant for Va ^ I'b mixing, such as that generated by the 
coefficients {g^^'^C)ab and {H^^C)ab in h^s- An analytic 
approximation to the electron number density inside the 
Sun is given by [52] n^/NA = 245e-i" '"^4i?,/7?,a ^ jg ^^g^f^l 

to define Ug = Ue — ^n„, a combination that often ap- 
pears in sterile-neutrino searches. This number density 
has a similar approximation, Hs/Na = 223e~^°'^*^/^®. 
The two linearly independent combinations can therefore 
be taken as GpUe ~ 1.32 x io-2Og-io.54ii/ii0 Q^y 
GfUs 0^ 1.20 X io-2Oe-io.54fl/i{0 Q^y^ corresponding 
to a neutron contribution of GpUn = 2Gp{ne — Ug) — 
0.24 X io-20g-io.54i?/_Ro ^j^g effective hamilto- 

nian. These quantities set the scale for matter effects in 
the Sun. 

Next, consider the detection of atmospheric neutri- 
nos. Upward-going neutrinos pass through the Earth 
and therefore experience higher matter potentials than 
the downward-going neutrinos, which traverse the less 
dense atmosphere and a small amount of bedrock on 
their way to the detector. A crude estimate of the mat- 
ter potential in this case can be obtained by assum- 
ing that the Earth consists of roughly equal numbers 
of protons, neutrons, and electrons. Using the aver- 
age number density then yields the approximate value 
GpUe — GpUn — 1.5 X 10~^^ GeV. This produces a mat- 
ter potential similar to that from the Sun at R/Rq ^ 2/5. 

Overall, the contribution to /leff from matter ranges 
from about lO^^o GeV to IQ-^s GeV. This means that 
matter effects must be incorporated when the contri- 
butions from mass or Lorentz violation lie near these 
values. This range is comparable to the scale of coef- 
ficients for Lorentz violation that originate as suppressed 
effects from the Planck scale. Note also that most terres- 
trial experiments involve neutrinos that traverse at least 
some amount of bedrock or other shielding materials, 



which can result in substantially different conventional 
or Lorentz-violating dynamics relative to the vacuum- 
oscillation case [53]. 



D. Neutrino oscillations 

The analysis of neutrino mixing proceeds along the 
usual lines. The effective hamiltonian can be diagonal- 
ized with a 6 x 6 unitary matrix UeS'- 

heS = UlgE,sU,B, (19) 

where Ees is a 6 x 6 diagonal matrix. In contrast to the 
Lorentz-covariant case, where mixing without sterile neu- 
trinos involves only three propagating states, here mix- 
ing without sterile neutrinos may occur with six states. 
This means that there can be up to five energy-dependent 
eigenvalue differences for Lorentz-violating mixing, re- 
sulting in five independent oscillation lengths instead of 
the usual two. 

Denoting the six propagation states by the amplitudes 
Bj{t;p) with J = 1,...,6, we can write Bj{t;p) = 
Ujaba{t;p) + Ujada{t;p), where we have split UeS into 
6x3 matrices UcB = {U ,U). The time evolution opera- 
tor may then be written as 

Sab{t) = {UlBe-'''^'''U,s)ab 

J \ <J jaUjb U JaUjb J 

where Ei^j) are the diagonal values of EcS- 

The probabilities for a neutrino of type h oscillat- 
ing into a neutrino or antineutrino of type a in time 
t are therefore Pi,^^i,^{t) = \S^^^^{t)\'^ or P^^^p„(t) = 
\Sii P, respectively. Similarly, for antineutrinos we 
have P„,^,,At) = \S.M.t}\^ or P,^^,Jt) = \S,^,M''- 
In terms of the matrices U and U, the probabilities are 



P,,^,„ (t) = 5„6 - 4 ^ Re {U*jJJjbUKaU*Kb) sin' + 2 ^ Im {U}„UjbUKaU*Kb) sin AjKt , (21a) 











Pu,^^M 




►.a(i) 



2 

J>K .]>K 



. 2 AjKt 



2 ^ Im {U*jJJji,UKaUKb) sin Ay^t , (21b) 



J>K J>K 



-4 ^ Re {U*j,UjbUKaU*Kb) sin' + 2 ^ Im {U*j,UjbU KaU*Kb) sin ^jkI , (21c) 



2 

J>K J>K 



. 2 AjKt 
mn 

2 



+ 2 ^ Im {U*jJJjbUKaUKi,) sin AjKt , (21d) 



J>K J>K 
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where the effective-energy difference is denoted by 

E. CPT properties 

With a conveniently chosen phase, CPT may be im- 
plemented by the transformation 

(22) 

This yields precisely the expected result when applied to 
hcff: the CPT-conjugatc hamiltonian h^g'^ = cr'^hlga^ 
can be obtained from Eq. (14) by changing the sign of 
the CPT-odd oj, and g coefficients. Then, h^g^ = hes 
when aL and g vanish, as expected. A notable feature 
here is that independent mass matrices for neutrinos and 
antineutrinos cannot be generated as has been proposed 
[54]. Greenberg has recently proved that this result is 
general [5]. 

Under CPT, the transition amplitudes transform as 
S...,{t) ^ Sl,^{-t), (23a) 

S...dt) ^ -S:^,A-t). (23b) 

These relations become equalities if CPT holds. The first 
relation then yields the usual result, 

CPT invariance =^ P^,^^„ (i) = Pp^^p,(t). (24a) 

This property has long been understood and has been 
identified as a potential test of CPT invariance [25]. How- 
ever, the negation of terms in this result produces a state- 
ment that may be false in general because CPT viola- 
tion need not imply Py^^y^(t) ^ P^^^vb{t)- Examples 
of models that violate CPT but nonetheless satisfy Eq. 
(24a) are given in Sec. IV. 

The above property addresses the relationship between 
v ^ v and v ^ D mixing. There is also an analogous 
property associated with ;/ i/ mixing. Thus, for CPT 
invariance, relation (23b) yields the additional result: 

CPT invariance =^ P^^^^^{t) = P^^^^^{t). (24b) 

This property may also provide opportiuiitics to test for 
Lorentz and CPT invariance. Note, however, that nega- 
tion of its terms produces a statement that may be false 
in general, as in the previous case. 

Finally, we emphasize that the presence of CPT vi- 
olation increases the number of independent oscillation 
lengths without the addition of sterile neutrinos. In the 
general case, nonzero coefficients for CPT violation in the 
effective hamiltonian (14) can generate up to six indepen- 
dent propagating states, rather than the usual three. 

F. Reference frames 

The presence of Lorentz violation makes it necessary 
to specify the frame in which experimental results are 



reported. Coordinate invariance of the physics, in par- 
ticular observer Lorentz invariance [2], ensures that the 
analysis and measurements of an experiment can be per- 
formed in any frame of reference. However, it is conve- 
nient to have a standard set of frames to facilitate com- 
parisons of different experiments. In the literature, mea- 
surements are conventionally expressed in terms of co- 
efficients for Lorentz violation defined in a Sun-centered 
celestial equatorial frame with coordinates (T, X, Y, Z) 
[55]. For our present purposes, it suffices to identify the 
Z direction as lying along the Earth's rotational axis and 
the X direction as pointing towards the vernal equinox. 
The coefficients for Lorentz violation in any other iner- 
tial frame can be related to the standard set in the Sun- 
centered frame by an observer Lorentz transformation. In 
general, this transformation includes both rotations and 
boosts, but boost effects are frequently neglected because 
they introduce only terms suppressed by the velocity (3 
between frames, which is typically < 10~^. Recently, 
studies of some /3-suppressed terms have been performed 
in the context of high-precision clock-comparison exper- 
iments [15, 16] and resonant cavities [21, 22]. 

The existence of orientation-dependent effects makes 
it useful to define a standard parametrization for the di- 
rection of neutrino propagation p and the corresponding 
ei, £2 vectors in the Sun-centered frame. A suitable set 
of unit vectors is given by 

p = (sin6cos<l?,sin6sin<l>, COS0), 
ei = (cos©cos$,cos9sin$, — sin6), 
€2 = (-sin$,cos$,0), (25) 

where & and $ are the celestial colatitude and longitude 
of propagation, respectively. We remark that these quan- 
tities are related to the right ascension r and declination 
d of the source as viewed from the detector by = 90°-|-rf 
and $ = 180° -t- r. 

In the remainder of this subsection, we provide some 
technical comments about the frame-dependence of our 
choice of spinor basis in Sec. II B. This basis is nor- 
mally associated with massless fermions, so the pres- 
ence of mass or Lorentz violation means that even with 
a covariant normalization the corresponding amplitudes 
are no longer scalar functions under observer Lorentz 
transformations and hence are frame dependent. How- 
ever, our basis suffices for perturbative calculations in 
which the physically significant states are affected only 
by masses and coefficients for Lorentz violation that are 
small relative to \p\ , while the complexity of the general 
Lorentz-violating case makes the decomposition into a 
covariant basis impractical. Moreover, despite the frame- 
dependent nature of the calculation, the probabilities 
(21) are frame independent at leading order. In the usual 
case, frame independence follows from the Lorentz-vector 
nature of the exact 4-momenta {E(^jy,p), which implies 
the products E(^j)t—p-x are Lorentz scalars, and from the 
constancy and frame-independence of the mixing matrix 
Ucs. It turns out that a version of these properties holds 
in the present case, as we show next. 
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First, we observe that the elements of the 6x6 matrix 
If'K^cff ^ \p\) are scalars under observer Lorentz trans- 
formations at leading order in small quantities. Next, 
note that the matrix Ues diagonalizes |p|(/ieff — so 
its elements can be chosen to be observer Lorentz scalars 
as well. In turn, this means that the diagonal elements 
\p\ {E(^j) — are also observer Lorentz scalars, since they 
are functions of the elements of |pl(/ieff ~ IpI)- Prom this 
result, it follows explicitly that the neutrino dispersion 
relations E'^j^ — p'^ are observer Lorentz scalars at lead- 
ing order, since 

Elj)-f = (E^j) + \mE(j)-\p\) 

2\p\iE^j^-\p\). (26) 

The 4-momentum is therefore a vector under observer 
Lorentz transformations to leading order, as desired. 
Combining this property with the scalar character of Ucff 
implies that the leading-order transition amplitudes and 
probabilities (21) are covariant under observer Lorentz 
transformations, as claimed. 



III. SENSITIVITIES 
A. Existing constraints 

To date, there is no compelling experimental evidence 
for nonzero coefficients for Lorentz violation in any sec- 
tor. Theoretical predictions of the size of the effects de- 
pend on the underlying model. However, the natural 
scale for a fundamental theory is the Planck mass mp, 
which is about 17 orders of magnitude greater than the 
electroweak scale mw relevant to the SM and roughly 
30 orders of magnitude greater than the scale of neu- 
trino masses, if they exist. It is plausible that any 
observable Lorentz- and CPT-violating effects are sup- 
pressed by one or more powers of the dimensionless ratio 
r = m/mp < 10~^^, where m is the relevant low-energy 
scale and top is the Planck mass [7]. In contrast, the 
scale of observed neutrino oscillations is < 0.1 eV, which 
enters as a squared mass Am^ < 10~^° GeV. At physi- 
cally relevant energies, lO"'^ GeV < E < 10^ GeV, the 
oscillation physics is determined by the dimensionless ra- 
tio Tlx = Am'^/E'^. Remarkably, the two dimensionless 
ratios r and r^, have a similar range, so the natural size 
of Lorentz- and CPT-violating effects may be comparable 
to the natural size of neutrino-oscillation effects. 

Certain experiments in the fermion and photon sectors 
have achieved sensitivities corresponding to dimension- 
less suppressions of roughly 10~^°. Since the coefficients 
for Lorentz violation in the various sectors can be related 
either directly through symmetries or indirectly through 
radiative corrections, it might seem that existing exper- 
imental constraints severely restrict the possibilities for 
Lorentz violation in neutrinos. In fact, this expectation 
is incorrect, as we discuss next. 



In the context of hcs, the relevant coefficients are 
(^L)ab arid {cl)'^'i^ , since these appear directly in the 
charged-fermion sector of the SME. A decomposition of 
the multi-flavor QED limit of the charged-lepton sector 
can be performed in analogy with Eq. (2). It produces 
the identification 

Kb = |K-«h)^6. 

C = ^(0,-0^):^, (27) 

where (cp)^^ and (ap)J^(, are coefficients in the SME that 
couple to right-handed Icptons and therefore leave unaf- 
fected the active neutrinos at tree level. On this basis, it 
might naively appear that the charged sector is sensitive 
to more combinations of coefficients for Lorentz violation 
than the neutrino sector. However, the mass hierarchy of 
the charged leptons e, /U, r implies that only coefficients 
that are diagonal in flavor space appear in leading-order 
perturbative calculations. As a result, e, fi, t effectively 
decouple, resulting in three independent copies of the 
fermion sector in the Lorentz- and CPT-violating QED 
extension. This implies that unsuppressed sensitivity to 
Lorentz violation in the charged-lepton sector involves 
only flavor-diagonal components. Moreover, the decou- 
pling also implies that certain coefficients such as a^^, 
a^^, o^^ are physically unobservable, further reducing 
the total number of coefficients affecting charged leptons. 
Taken togethc^r, these factors ensure that the CPT-odd 
sectors of charged leptons and neutrinos are completely 
independent at tree level. Similar arguments apply to 
parts of the CPT-even sector as well. We therefore con- 
clude that neutrinos are sensitive to a greater number of 
coefficients for Lorentz violation than the charged lep- 
tons, and at tree level most of these coefficients are inde- 
pendent from those accessible with e, or r leptons. 

Particularly stringent constraints exist on some com- 
ponents of the charged-lepton coefficients and 
Although these arc linearly independent of neutrino- 
sector coefficients at tree level, it is natural to ask 
whether radiative corrections to these components can 
be used to constrain possible neutrino effects. As an 
example, Ref. [28] explores the possibility that eV-size 
effects in heavy sterile neutrinos could evade the con- 
straints in the charged-lepton sector, finding that within 
a standard seesaw mechanism the existence of large h^- 
type coefficients for sterile neutrinos tends to produce 6^ 
coefficients in the charged-lepton sector that conflict with 
observation. In this work, we neglect seesaw-induced co- 
efficients because they are suppressed by the large-mass 
scale. However, it is of interest to ask whether radia- 
tive corrections alter the tree-level independence of the 
charged- and neutral-lepton sectors. 

For simplicity, we restrict our discussion to the relevant 
0** and coefficients, although related remarks apply 
also to c^^ and d^^^ coefficients. The leading-order radia- 
tive corrections are linear in the coefficients for Lorentz 
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violation. However, loops involving weak-interactions arc 
heavily suppressed by additional factors at the relevant 
energies, while strong interactions play no role. We can 
therefore restrict attention to the QED extension. In this 
case, general properties of the coefficients for Lorentz vi- 
olation under the discrete symmetries C, P, and T im- 
ply that corrections to coefficients involve only other 
6^ type coefficients [56]. As a result, although the con- 
straints from chargcd-lepton experiments may restrict &^ 
in the neutrino sector of the SME, the a'^ coefficients 
are unaffected and so oj, is unconstrained. Thus, the 
independence of the charged- and neutral-lepton sectors 
remains valid for radiative corrections. 



B. General features 

In the presence of Lorentz and CPT violation, a wide 
range of unconventional neutrino behaviors can occur. 
These include unusual energy dependence, direction- 
dependent effects, and neutrino- antineutrino mixing. 
Specific examples of these behaviors are illustrated in the 
examples presented in Sec. IV. Here, we focus on some 
general features of experimental sensitivities to Lorcntz- 
and CPT-violating effects. Some of these have been dis- 
cussed in the context of the minimal SME in our earlier 
work [29], but the present discussion holds for the full 
theory (14) and generically for operators of nonrenor- 
malizablc dimension. 

Figure 1 shows an estimate of the coverage in base- 
line distance L versus energy E of the currently pub- 
lished neutrino-oscillation data. Included in the evidence 
for oscillations are observations of solar neutrinos by Gl- 
and Ga-based experiments [37-40], Super Kamiokande 
(SK) [41], and SNO [42]; and of atmospheric neutrinos by 
SK [36] , reactor-based KamLAND [43] , and accelerator- 
based LSND [44] and K2K [45]. Null results include the 
reactor experiments Bugey [57], GHOOZ [58], Gosgen 
[59] , Palo Verde [60] , and various accelerator-based short- 
baseline experiments including, for example, the high- 
energy experiments BNL-E776 [61], GGFR [62], GHO- 
RUS [63], NOMAD [64, 65], NuTeV [66], and the low- 
energy KARMEN [67]. A number of new accelerator- 
based experiments are likely to produce interesting re- 
sults in the near future. These include the short-baseline 
(L ~ 500 m, S ~ 1 GeV) MiniBooNE experiment [68] de- 
signed to test the LSND anomaly, and the long-baseline 
{L ~ 700 km, £; ~ 1 GeV) ICARUS [69], MINOS [70], 
and OPERA [71] experiments, which are planned to test 
the atmospheric-oscillation hypothesis. Also shown on 
the figure are the approximate effective regions associated 
with the matter potentials for the Sun and the Earth. 

The unusual energy dependence can be viewed as 
a consequence of the dimensionality of the coefficients 
for Lorentz violation. The standard scenario for neu- 
trino oscillations involves mass-squared differences Am^ 
that combine with the baseline distance L and the 
neutrino energy E to yield the physically relevant di- 




E (GeV) 

FIG. 1: Approximate sensitivities of various experiments. 
Lines of constant L/E (solid), L (dashed), and LE (dot- 
ted) are shown, giving approximate sensitivities to the quan- 
tities {m,m5}, {a", b'',H and {c'"' , d"" , g""" } , respec- 
tively. Also shown are the approximate effective regions for 
the matter potential in the Sun and Earth. 

mensionless combination Am^L/E. However, Eq. (14) 
shows that Lorentz-violating oscillations generated by 
the dimension-one coefficients a^, b^, H^" are controlled 
by the dimensionless combinations a^L, b^L, H^'^L, 
while those generated by c'^" , d'^" , g^'^'^ are controlled 
by c^'^LE, d^'^LE, g^""^LE. More generally, oscillations 
generated by a coefficient (kd)^'" for a Lorentz-violating 
operator of nonrenormalizable dimension n = d + 5 are 
controlled by {kd)^-LE'^. 

Figure 1 illustrates these various energy dependences. 
Lines of constant L/E, L, and LE are plotted, bounding 
approximate regions of experimental sensitivity to con- 
ventional mass-squared differences, dimension-one coef- 
ficients, and dimensionless coefficients, respectively. For 
each nonzero coefficient in h^a, a bounding line on this 
figure exists above which the corresponding Lorentz- 
violating effects become of order one. Given such a line, 
any experiments locatc;d near or above it can be affected 
by the associated coefficient, but experiments below it 
have limited or no sensitivity. For example, the region of 
limiting sensitivity for a hypothetical dimensionless coef- 
ficient of magnitude ~ 10~^* is bounded approximately 
by the dimensionless line satisfying LE ~ 10^*, which is 
the dotted line running just below KamLAND. Exper- 
iments lying above this line, such as KamLAND, SNO, 
and SK, coidd be sensitive to the effects of this coefficient. 
Note that approximate regions of experimental sensitiv- 
ity to coefficients (fc^)^' ' of dimension 1 — d could also be 
identified on the figure. They would be bounded by lines 
of constant LE''' with d> 1, which have negative- integer 
slopes. 

Figure 1 also reveals that experiments and data al- 
low probes well below the 10^^^ Planck-suppression level. 
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For instance, the various null results from short-baseline 
reactor and accelerator experiments could be reanalysed 
to yield upper bounds on certain coefRcients for Lorentz 
violation. Thus, the high-energy experiments CHORUS 
and NOMAD found no evidence of i'e,ij, at energies 

E ~ 100 GeV and at distances L ~ io^^ GeV"\ which 
suggests that reanalyses of these experiments would yield 
interesting new sensitivities of roughly 10~^^ GeV to 
dimension-one coefficients and roughly 10"^*^ to dimen- 
sionless coefficients. A similar situation holds for low- 
energy experiments such as CHOOZ, Palo Verde, and 
KARMEN in the Ve sector. From Fig. 1 we see that, 
relative to CHORUS and NOMAD, CHOOZ and Palo 
Verde might be expected to have comparable sensitivi- 
ties to dimension-one coefficients but reduced sensitiv- 
ity to dimensionless ones, while KARMEN has compa- 
rable sensitivity to conventional mass e£Fc!cts. In each 
case, the attainable sensitivities also depend on various 
experiment-dependent factors, so individual reanalyses 
are required to make definitive statements. 

Another unusual effect due to Lorentz violation is 
direction-dependent neutrino behavior, a consequence of 
rotation-symmetry violation. This has consequences for 
comparisons of results between different terrestrial exper- 
iments or for the analysis of experiments involving mul- 
tiple sources, since the orientation of the neutrino beam 
or the location of the source relative to the detector can 
affect neutrino oscillations. Rotation-symmetry violation 
also implies that the daily rotation of the Earth about its 
axis induces apparent periodic changes of the coefficients 
for Lorentz violation in the laboratory, which would be 
manifest as temporal variations in neutrino oscillations. 
These variations occur at multiples of the sidereal fre- 
quency io^ (23 h 56 min). Similarly, in the presence 
of rotation-symmetry violation, neutrinos emitted from 
the Sun in different directions undergo different oscilla- 
tions, which may produce observable annual variations 
arising from the change in the location of the detector 
as the Earth orbits the Sun. All these temporal varia- 
tions with appropriate periodicity provide unique signals 
of Lorentz violation in neutrino oscillations. Moreover, 
they can also yield interesting sensitivities to certain co- 
efficients. For instance, SK found that the shape of the 
solar-neutrino flux matches the expected value to within 
about 5% over the year [41]. The Sun-Earth distance 
is L ~ 1027 GeV-i,' and LE ~ 10^5 for the SK energy 
range, so a reanalysis of the SK data might achieve im- 
pressive sensitivities of ^ 10"^^ GeV to dimension-one 
coefficients and ~ 10"^^ to dimensionless ones, compa- 
rable to the best experimental sensitivities achieved for 
other sectors of the SME. 

Another interesting feature of Lorentz violation in- 
volves novel resonance effects in neutrino oscillations. In 

the conventional case with neutrino masses, the usual 
MSW resonances [35] arise when the local matter envi- 
ronment is such that neutrino interactions become com- 
parable to mass effects, thereby drastically changing the 
character of the hamiltonian. The presence of Lorentz 



violation can trigger several other types of effects, in- 
cluding resonances without mass or matter that involve 
different coefficients for Lorentz violation, resonances in- 
volving coefffcients for Lorentz violation and mass terms, 
resonances involving coefffcients for Lorentz violation and 
matter effects, and various combinations of the above. 
The earliest example of an explicit vacuum resonance 
in a two-generation model involving a mass term and 
a single nonzero coefficient (ol)'^ for Lorentz and CPT 
violation is given in Rcf. [25] . An example of a vacuum 
resonance in a three-generation model involving two coef- 
ficients (ol)^ and (cl)"^"^ for Lorentz and CPT violation 
occurs in the bicycle model of Ref. [29] . We emphasize 
that resonances due to Lorentz violation can occur in the 
vacuum as well as in matter, and not only at particular 
energies but also for particular directions of propagation. 
Note also that, even away from the resonance regions, 
matter effects may be important when considering mass 
terms or coefficients for Lorentz violation that have lines 
of sensitivity near or above the Sun- or Earth-potential 
regions shown in Figure 1. 



C. The LSND anomaly 

In the LSND experiment [44] , copious numbers of neu- 
trinos were produced from the decay of 7r+ at rest. This 
process is dominated by the decay 7r+ A^'''^'^ followed 
by /i^ e^VeV^. A small excess in Ve was seen, inter- 
preted as the oscillation with a small probabil- 
ity of about 0.26%. This result is difficult to accommo- 
date within the context of the conventional global anal- 
ysis [31], in which two mass-squared differences are used 
to describe solar and atmospheric oscillation data. The 
solar data appear consistent with a mass-squared differ- 
ence 6vp? ~ 10~^ eV^, while the atmospheric data sug- 
gest a second mass-squared difference Am^ 10~^ eV^. 
The regions of limiting sensitivity to these mass-squared 
differences are shown in Fig. 1, where lines of constant 
LjE with values LjE ^ W'' GcV^^ and LjE ^ lO^i 
GeV~2 can be seen. Experiments lying significantly be- 
low these lines, including LSND, should be insensitive 
to oscillations caused by ?>w? and Am^. This illustrates 
the difficulty in explaining the LSND result within the 
conventional framework without introducing additional 
mass-squared differences. 

A resolution of this LSND anomaly without the intro- 
duction of sterile neutrinos might emerge from the un- 
usual energy dependence, the directional dependence, or 
the neutrino-antineutrino mixing introduced by Lorentz 
violation. For example, equal numbers of i^^, j/g, and v^, 
are produced in LSND, so if i/g mix with v,, then the ob- 
served excess in i/g may be a result of z/g <-> Ve mixing 
rather than ^ mixing. We note, however, that if 
the possible direction dependence is neglected then Fig. 
1 shows that a simple solution based either on the un- 
usual energy dependence ox an v v mixing is likely 
to be hindered by existing null results in the sector. 
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from low-cncrgy experiments such as CHOOZ and Palo 
Verde or from high-energy experiments sueh as CHO- 
RUS, NOMAD, and NuTeV. Indeed, from this figure we 
see generically that to explain the LSND result one needs 
a mass-squared difference of about 10~^^ GeV^ = 10~^ 
eV^, a dimension-one coefficient of about 10^^* GeV or 
a dimensionless coefficient of about 10~^^. Note that 
each of these has consequences for other experiments, 
depending on flavor content. For example, the upcom- 
ing MiniBooNE experiment is designed to test the same 
oscillation channel and will therefore be sensitive to all 
three possibilities. 

IV. ILLUSTRATIVE MODELS 

To illustrate some of the novel behaviors of neutrino 
oscillations in the presence of Lorentz and CPT viola- 
tion, we next consider a number of simple special cases 
of the theory (14) with only one or a few nonzero co- 
efficients. For each case, some of the ways that the 
unusual neutrino behaviors might affect current obser- 
vations are quantitatively examined. Also, we simplify 
expressions by adopting temporary notation for the spe- 
cific nonzero coefficients for Lorentz violation within each 
model: quantities carrying a ring accent, such as c, de- 
note rotation-symmetric coefficients; while those with a 
hacek accent, such as c, denote anisotropic coefficients. 

A. Rotationally invariant models 

The rotation-invariant restriction provides an interest- 
ing special limit of the theory (14). While difficult to 
motivate without knowledge of the underlying mecha- 
nism leading to Lorentz and CPT violation, rotation- 
invariant or so-called 'fried-chicken' (FC) models are at- 
tractive because rotation symmetry can significantly re- 
duce the complexity of calculations, thereby providing a 
simple context within which to study the unusual neu- 
trino behaviors arising from Lorentz violation. 

Restricting hes to FC terms leaves only four matri- 
ces, {m^)at, {a^t,, {cl)Z and {clYJI = ^('-^lYL^"'- As 
described in Sec. II B, the trace (cl)°^ — {clYJi) unob- 
servable and may be set to zero, so only three of these 
matrices are independent. Dropping the irrelevant ki- 
netic term and assuming rotation invariance in the Sun- 
centered (T, X, Y, Z) frame for definiteness, the 6 x 6 ef- 
fective hamiltonian reduces to the block-diagonal form 

(M^f = di^g[{my {2E) + {a^f -^{clY^E)^^, 
(mV(2i^)-K)^-|(c^)^^i?):j.(28) 

This hamiltonian provides a general FC model of three 
active neutrinos. The generalization to additional light 
or massless sterile neutrinos is straightforward. 

With the exception of the original proposal for Lorentz 
violation in neutrinos [2] and the recent work in Ref. [29] , 
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FIG. 2: Lines of limiting sensitivity for 5m^ ~ 10 ^ eV^, 
Am^ ~ 10"^ eV^ 5c ~ 10-2^ Ac ~ 10"". The shaded 
regions are those of Fig. 1. 



which address both rotation-invariant and anisotropic ef- 
fects with and without CPT violation, existing works on 
the subject [24-27] involve limited special cases of the 
general FC model (28). The bulk of the literature re- 
stricts attention to the two-generation special case and 
neglects either the (a/,)-^ term or the (cl)^^ term. A 
plethora of unexplored models and effects exists. 

It might seem logical to impose spherical symmetry in 
a special frame such as the cosmic microwave background 
(CMB) frame. However, if rotation symmetry is assumed 
in the CMB frame then the coefficients in Eq. (28) differ 
from (cl)-^, (cl)-^-^ in the standard Sun-centered frame, 
being instead (ol)", (cl)°° in the CMB frame. Relat- 
ing the latter to the standard Sun-centered frame or any 
other experimentally attainable frame introduces direc- 
tion dependence due to the motion of our solar system 
in the CMB frame. The relevant hamiltonian then also 
involves spatial components of the coefficients, so it dif- 
fers from Eq. (28) and is instead an anisotropic limit of 
the theory (14). 

Although the FC model (28) is rather limited consid- 
ering the wealth of possible effects contained in the full 
theory (14), and although it has little theoretical motiva- 
tion other than calculational convenience, further study 
of this model is useful because it provides a readily work- 
able context within which to gain insight about possible 
signals of Lorentz and CPT violation. This is illustrated 
in the few simple examples discussed in this subsection. 



1. Example: {cLjab 

A particularly simple FC model consists of a single 

nonzero coefficient matrix such as (c/,)^^. Some features 
of this model are similar to the conventional massive- 
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neutrino case, but there is unusual energy dependence. 
Here, wc take advantage of this energy dependence to 
illustrate one type of mechanism through which Lorentz 
violation might provide a solution to the LSND anomaly. 

Lines of limiting sensitivity for the two mass-squared 
differences Sm^ and Ato^ used in the conventional global 
analysis are shown in Fig. 2. The mixing angles are such 
that fe oscillations are almost completely controlled by 
Sm^. Therefore, one can expect to see only z/g mixing in 
KamLAND, in solar neutrino experiments, and possibly 
in the lowest-energy atmospheric-neutrino experiments. 
The observed atmospheric oscillations are due to Am?, 
which controls U/j, ^ Ur mixing. Since LSND lies well 
below both the 5m^ and the Am^ lines, no oscillations 
are predicted. 

Replacing the mass-squared differences dm'^ and Am^ 
with a nonzero coefficient matrix {cl)^'^ produces an ef- 
fective hamiltonian h^s that can be parametrized as de- 
scribed in Appendix B, using two eigenvalue differences 
and CKM-like mixing angles and phases. For simplic- 
ity, we choose here to mimic the usual solution by taking 
vanishing phases and O^^ , and we consider only the case 
— T^l^- This leaves three degrees of freedom: two 
eigenvalue differences, and one mixing angle = ^■ 
It turns out to be convenient to work with two linear 
combinations of the eigenvalue differences, defined by 

5c ^ |((c.)f3y-(c,)f,^), 

Ac ^ |((c^)^-(c^)fiP. (29) 

The probabilities for this case are then 

- s\r? e sin2((Ac + 5c^LE/2) 

-cos^esm'^{5cLE/2), 
P,^^,^ = P,^^,^^^sm'^2esm'^{AcLE/2), 
P.,^u^ = -\sm^2esm^{AcLE/2) 

+ sin^ e sin2((Ac -|- 5c)LE/2) 

+ cos^esm^{6cLE/2). (30) 

The corresponding antincutrino expressions are identical. 

A possible approach is illustrated in the figure. The 
line of sensitivity for the larger difference Ac can be cho- 
sen to lie just above CHOOZ and LSND. This produces 
only a small effect in these experiments and may pro- 
vide an explanation for LSND that may not conflict with 
CHOOZ. The remaining difference 5c can then be chosen 
to explain atmospheric data. The above situation some- 
what resembles the conventional mass solution, with the 
role of 5m?'/2E replaced by AcE and that of Am? /2E 
replaced by 5cE. The angle 6 is the analogue of the solar- 
neutrino mixing angle. However, the energy dependences 
of the two cases differ substantially, as is also evident 
from the figure. 

To explore quantitatively how this approach might 
work, consider the positive LSND and KamLAND re- 



sults. KamLAND detects v,, from distant reactors and 
found about a 61% reduction in the flux. Most reac- 
tors are 138-214 km from the detector, and the corre- 
sponding Vf. energies fall in the range 1 MeV < E < 
10 MeV. If KamLAND lies well above the Ac line, 
the relevant quantity is the average survival probabil- 
ity (Pp^^p^) = 1 — ^ sin^ 20 ~ 61%, yielding a mixing 
angle given by sin^ 29 ~ 0.78. Also, assuming LSND 
is in a region of small oscillation effects, then we can 
approximate Pp^^p, w ^sin^ 20{AcLE/2)^ ~ 0.26%. 
Then, for E ~ 45 MeV and L ~ 30 m we obtain 
Ac ~ 2.4 X 10"^^. Thus, in this simple scenario, these two 
experiments suggest coefficient values near sin^ 29 ~ 0.78 
and Ac ~ 2.4 x 10"^'', in agreement with the estimates 
of Sec. inc. 

The remaining coefficient 5c can then be chosen to 
match observed atmospheric-neutrino effects. The coef- 
ficient Ac is relatively large in this region and generates 
rapid oscillations. Averaging over these for any value 
of 5c leaves a muon-neutrino survival probability of ei- 
ther P^^^^^ ~ 0.54 0.27sin2(,5cL£;/2) or P,.^^^,, ~ 
0.77- 0.73sin^((5cLi;/2), depending on the solution for 
9. Note that the latter expression resembles the usual 
maximal-mixing solution within an overall scale factor, 
except for the unusual energy dependence in the oscilla- 
tion length. 

Interestingly, atmospheric electron-neutrino oscilla- 
tions are present in this model but are largely unobserved 
due to a compensation mechanism. The averaged i^e sur- 
vival probability is Pi^^^i,^ = 61%, as above, and the 
t^e i^fi mixing probability is P^^^^^ = 19.5%. The ob- 
served flux of atmospheric electron neutrinos is a combi- 
nation of the survival flux and the appearance flux from 
mixing with muon neutrinos. Since the ratio of muon 
neutrinos to electron ncTitrinos is approximately 2, the 
predicted effective flux of atmospheric electron neutri- 
nos is approximately 61% + 2(19.5%) ~ 100% of the flux 
in the absence of oscillations, in agreement with indica- 
tions from existing data. Essentially, this compensation 
mechanism works because the disappearance probabil- 
ity 1 — P^^^v^ of electron neutrinos given by Eq. (30) is 
a factor of two greater than the appearance probability 
Pvc^v,i of muon neutrinos from mixing, resulting in no 
net suppression in the total observed electron-neutrino 
flux. 

The compensation mechanism per se is independent 
of Lorentz violation and can be applied whenever 1 — 
Pj^^_,j^^ « 2P^^^^^, including in the conventional mas- 
sive case. Note, however, that Monte Carlo calculations 
suggest the flux ratio increases dramatically above 2 for 
energies over about 10 GeV [72], so the compensation 
mechanism is likely to fail at higher energies. Note also 
that, in the case of the above Lorentz-violating model, 
the rapid oscillations at high energies also help to mask i/e 
oscillations. Although these rapid oscillations can change 
the overall flux, they also tend to smooth away the ob- 
servable E and L dependences that form the basis for 
some analyses. 
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This simple model serves to illustrate a possible strat- 
egy that might remedy the eonflict between LSND and 
reactor experiments, but it may well introduce other con- 
flicts between LSND and accelerator experiments testing 
Vf. Vr and t'^ Vr [63, 64] or i'^ Vf, [65, 66]. 
Note also that some work has been done to check for 
unconventional energy dependences in the atmospheric 
data [73], suggesting that the usual energy dependence 
is preferred. However, these analyses are limited to two 
generations and do not consider possible direction depen- 
dences or mixing. A complete treatment would 
also need to include the effects of the Earth's matter po- 
tential, which introduces additional energy dependence. 
The point is that Gpn,, ^ 10~^^ GeV for the Earth, 
and at atmospheric-neutrino energies this is comparable 
to the contribution from 8c shown in Fig. 2. In any 
case, interesting sensitivities to Lorentz violation could 
be achieved with a complete analysis of existing data. 



2. Example: (ai^f^^ ^ 0, (ci)JJ ^ 



We turn next to an FC model with mixed energy 
dependence, incorporating only two nonzero coefficients 
(a^)^^ = a and |(cl)JJ = c and no mass terms. This 
case includes both Lorentz and CPT violation but re- 
mains rotation symmetric. The presence of both a dimen- 
sionless coefficient and a dimension-one coefficient leads 
to unusual energy behavior in the vacuum-mixing angles 
as well as the oscillation lengths. This contrasts with the 
previous case, in which only the oscillation lengths have 
unconventional energy dependence. Note that both a and 
c are arbitrary to an unobservable phase, and therefore 
they can be taken real and nonnegative without loss of 
generality. 

The behavior in this model can be understood quali- 
tatively as follows. At low energies E a/c relative to 
the critical energy a/c, the a term dominates the effec- 
tive hamiltonian. As a result, Vr decouples from and 
so only t'g Vfj, mixing occurs. In contrast, for high 
energies E ^ a/c, c dominates and only t'^ ^ i^r mix- 
ing occurs. At intermediate energies E ~ a/c, the two 
terms are comparable and produce complicated energy 
dependence with mixing between all three neutrinos. 

This behavior is similar to the observed energy depen- 
dence in the solar-neutrino flux. In the usual analysis 
with massive neutrinos, the observed energy dependence 
is explained through matter effects. However, the same 
type of behavior can appear in Lorentz- violating scenar- 
ios even without matter. To demonstrate this, we need 



the probabilities for the current model: 

P^^^^^ = 1- 4:31-0^ 9 cos^esm'^inL/Lo) 

-sin'*6'sin2(27rL/Lo), (31a) 

P,^_,^ = 1 - sin2(27rL/Lo), (31b) 

P^^^^^ = l-4sin2 6'cos2 6lsin2(7rL/Lo) 

-cos^esm^{2TTL/Lo), (31c) 
= sin2^sin2(27rL/Lo), (31d) 

P^^„^^ = sin^^cos^6'(4sin^(7rL/io) 

-sin2(27rL/Lo)), (31e) 

P,^„,^ = cos2esin2(27rL/io), (31f) 

where 

sin2^ = a^/{a^ + c^E^), 
2n/Lo = ^/dFT^. (32) 

The antineutrino probabilities are again identical since 

the quantities sin^ 9 and Lq are symmetric under a — > 
—a. We remark in passing that this model serves as an 
example in which CPT is violated but the traditional test 
of CPT discussed in Sec. II E fails as an indicator of the 
CPT violation. 

The solar-neutrino vacuum-oscillation survival prob- 
ability is given by Eq. (31a). As usual, depending on 
the size of the coefficients, matter effects can drastically 
alter the survival rates. Consider, for example, a sim- 
ple matter-dominated case where the matter potential 
at the point of Ve production dominates heg. Assum- 
ing adiabatic propagation, neutrinos are produced in the 
highest-eigenvalue state of hcB{R ~ 0) and emerge from 
the Sun in the highest-eigenvalue state of hcff{R = Rq). 
The overlap between this state and an electron-neutrino 
state is proportional to sin^/\/2. Consequently, the av- 
erage survival probability for the matter-dominated case 
in an adiabatic approximation is 

(P^^^^J adiabatic = 5 sin^ 9. (33) 

In contrast, the average for the case where matter effects 
can be neglected is 

(P.e--e ) vacuum = 1-2 sin^ ^ + I sin^ 9. (34) 

These probabilities are plotted on Fig. 3 as a function of 
energy in units oi a/c. 

The observed flux is consistent with the figure, since 
low-energy experiments suggest an approximate survival 
probability of 1/2 [38-40], while higher-energy experi- 
ments favor about 1/3 [37, 41, 42]. Note that both cases 
shown in Fig. 3 yield an average survival probability of 
1/3 at E = h/\p2,c. By choosing the ratio a/c to coin- 
cide with the peak of the solar ®B spectrum (-Epoak — 6.4 
MeV), this simple massless Lorentz- and CPT- violating 
model can be made to reproduce the gross features of the 
observed solar-neutrino flux. This corresponds to impos- 
ing a/c ~ 9 MeV. 
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FIG. 3: Solar-neutrino survival probability assuming adia- 
batic propagation (solid), and average survival probability for 
vacuum oscillations (dashed) . 

The above discussion only depends on the ratio of coef- 
ficients. To get a sense of the size of coefficients required 
in a reaUstic case, we can consider what KamLAND im- 
phes for a and c. Taking a representative neutrino to 
have energy E — 5 MeV and baseline L = 200 km and 
assuming that it oscillates no more than once, the ratio 
a/coi 9 MeV and the survival probability Pu^^u^ ~ 61% 
can be used to extract approximate values a ~ 7 x 1G~^^ 
GeV and c ~ 8 x 10^^". The lines of sensitivity for these 
values on Fig. 1 are approximately L ~ 10^^ GeV"''^ and 
LE ^ 10^'', passing just above KamLAND and inter- 
secting in the solar-energy region, thereby producing the 
energy dependence seen in Fig. 3. 

3. Example: (m^)eM ^ 0, (01.)^^ ^ 

As a variation on the above model, wc next consider 
a special FC case with nonzero mass {'irP)efj. = rh^ and 
coefficient {aL)J^T — ^' fo'^ Lorentz and CPT violation. 
This model has many qualitative features of the previ- 
ous one. At small energies, the mass m controls mixing 
between fe and u^, while at large energies a' dominates 
and produces mixing between and Vt- 

The probabilities for this model arc given by Eqs. (31), 
(33), and (34), but with the definitions 

sm^e = m*/{m* + 4:a'^E^), 
27r/Lo = y/{my2E)-^ + a'^. (35) 

The analysis of this model parallels the previous case. 
Indeed, Fig. 3 also holds for the solar-neutrino proba- 
bilities in terms of m and of, using the scale shown on 
the top axis. Applying the same arguments as before 
yields the ratio rr? jb! ~ 18 MeV and candidate values 
m2 ~ 7 X 10-6 eV2 and a' ~ 4 x IQ-^^ GeV. 

A key difference between this case and the previous 
o-c model is the asymptotic behavior of the oscillation 
length. In the a-c case, Lq ^ 2n/{cE) at high ener- 
gies. In contrast, the oscillation length in the present m- 
a! model approaches a constant at high energies, Lq 
27r/a'. Consider the consequences for atmospheric neu- 
trinos. Note that in the high-energy limit of both cases. 



sin^ 9^0 and so Pp^^vc ^ 0, in agreement with ob- 
servation. However, the first model with c ~ 8 x 10^^*^ 
gives Lq ~ 2it/{cE) ^ (15 km GcV)/E, whereas the sec- 
ond model with a' ~ 4 x lO'^^ GeV yields Lq ~ 3100 km. 
These differ from the usual massive-neutrino explanation 
of the atmospheric data, which has Am? ~ 3 x 10"^ eV^ 
and results in Lq = ATrE/Am^ ~ 800£; km/GeV. 

We emphasize that both this special model and the 
previous one involve only two degrees of freedom, whereas 
the usual massive-neutrino solution requires two mass- 
squared differences and at least two mixing angles. In- 
cluding additional coefficients for Lorentz violation can 
only add flexibility to the analysis. For example, one 
might consider a combination of the two examples above, 
which would have four degrees of freedom. With addi- 
tional freedom, it seems likely that an appropriate sim- 
ple Lorentz-violating scenario could be constructed that 
would reproduce most oscillation data. This also sug- 
gests that existing data analyses appear insufficient to 
exclude many forms of Lorentz and CPT violation, or 
even to distinguish between oscillations due to mass and 
those due to Lorentz violation. 



B. Direction-dependent and <-> 1/ mixing models 

Lorentz violation naturally allows directional depen- 
dence in oscillation parameters through the violation of 
rotation invariance. An interesting subset of direction- 
dependent models are those involving v v mixing via 
nonzero g^'^'^ and H'^'^ coefficients in the theory (14). 
In the general case, nonzero mixing represents 

one way to generate as many as five distinct oscillation 
lengths without incorporating sterile neutrinos. How- 
ever, we limit attention in this subsection to a simple 
model that reveals some key features oi v v mixing. 
For illustrative purposes, it suffices to consider mixing 
in only one neutrino species, say Ve ^ i^e- This case 
may nonetheless have physical relevance, since it implies 
significant effects on reactor experiments and solar neu- 
trinos and might possibly also shed light on the LSND 
anomaly. 



1. General one-species model 

The restriction to the two-dimensional v^-Ve subspace 
radically simplifies the form of the effective hamiltonian 
(14). Since the coefficients (TO^)ee and {cLjee are real, 
they lead to terms proportional to the identity that have 
no effect on oscillatory behavior and can therefore be ig- 
nored. Moreover, Eq. (4) implies that {H^^C)ab is anti- 
symmetric in generation space, so {H^^^C)ee = H'^^c = 0. 
Therefore, the most general single-flavor theory without 
mass differences is given by a 2 x 2 effective hamiltonian 
containing only the coefficients (az,)^e ^^'^ {g'^"'^C)f,e = 
g^'^c for Lorentz violation. Note that both these terms 
are CPT odd. 
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For this general single-flavor model, the probabilities 
are identical in form to those of the usual two-generation 
mixing case: 

P -=1 — P =1 — P- 

= sm^2es\r?2-KL/Lo. (36) 

However, the mixing angle and oscillation length can 
have nontrivial 4-momentum dependence. They are 
given by the expressions 

sin^2^ = fl + ^^ife^ \\ (37) 

Note that these can also be written directly in terms of 
the neutrino-propagation angles and $ defined in Eq. 
(25). 

2. Example: g^^c 7^ 

As an explicit example, we consider a maximal-mixing 
special case of the general single-flavor model for which 
the only nonzero coefficient is = g. In terms of 

the propagation angles Q and the oscillation length is 
found to be 

27r/Xo = |£^sine,g|, (38) 

and the mixing angle is sin^ 20 =1. As in the previ- 
ous examples, this case has unconventional energy depen- 
dence, but unlike previous examples it includes neutrino- 
antineutrino mixing and also dependence on the direction 
of propagation through the propagation angle &. 

To illustrate the effects of the direction dependence, 
consider atmospheric neutrinos detected in the SK detec- 
tor. Neutrinos that enter the detector from the celestial 
north or south have sin 6 = and therefore do not oscil- 
late. In contrast, neutrinos propagating in the plane par- 
allel to the Earth's equatorial plane have sin 8 = 1 and 
experience maximal mixing [74]. Analyses of SK data 
often neglect the difference between Ug and Pg, so they 
may be insensitive to this effect because the total flux of 
electron neutrinos and antineutrinos is unchanged. How- 
ever, the same type of directional dependence can arise in 
more complicated scenarios with t'e I'/n i't mixing, 
and this could drastically affect the up-down asymmetry 
measurements of SK. 

As another example consider KamLAND, which de- 
tects neutrinos from several reactors at different loca- 
tions. The total flux (^totai(i') of can be written 

4>toME) = 5^ <i>j{E)P-,^^^^{E, Lj, Sj), (39) 

3 

where the (pj (E) are the fluxes from the individual reac- 
tors in the absence of oscillations, and &j are appropriate 



propagation angles determined by the relative positions 
of the reactors and the KamLAND detector. We can ap- 
proximate the positions of the reactors as being located 
in the plane tangent to the surface of the Earth at the 
location of the detector. It follows that neutrinos from 
reactors positioned directly north and south of the de- 
tector have Qj ~ 180° — x and Qj ~ Xi where x — 36° 
is the latitude of the detector. In contrast, neutrinos 
arriving from the cast or west have ~ 90°. This re- 
sults in an approximate allowed range for the Qj given by 
sin^ 6j > sin^ Xi implying that the Pg from every reac- 
tor experience some degree of oscillation on their way to 
the KamLAND detector. However, the net result differs 
from the flux in a comparable rotation-symmetric model 
with a dimcnsionlcss coefficient. 

For solar neutrinos, the allowed range for Q is given by 
sin^ O > cos^ r] ~ 0.85 because the Earth's orbital and 
equatorial planes differ by approximately rj = 23°. The 
true value of sin^ oscillates between sin^ 6 = 1 in the 
spring or fall and sin^ ~ 0.85 in the summer or win- 
ter. This simple model therefore predicts a semiannual 
variation in the solar-neutrino data. 

As suggested in Sec. IIIC, oscillations of Vg into Pg 
may provide an alternative approach to resolving the 
LSND anomaly. If the LSND result is reinterpreted as 
an oscillation of Vg into Pg, then the transition probabil- 
ity is likely to be comparable to the reported value of 
about 0.26% because roughly equal numbers of Vg and 
Pf^ are produced. Since mixing in this model is caused 
by the dimensionless coefficient g, a reasonable strategy 
here is similar to that adopted for the 6c-Ac model in 
Sec. IV A 1, where a dimensionless coefficient is chosen 
to have its line of sensitivity just above CHOOZ and 
LSND in Fig. 1. This causes a small oscillation in LSND 
but avoids the null constraints from reactor experiments. 
Taking the energy of a typical Vg to be about £^ = 35 
McV and the distance to be L = 30 m in LSND, and as- 
suming that the small transition probability is due to a 
small L/Lo,we can write P^^^i^^ = sin^ 29 sin^ 2'!rL /Lq ~ 
{smOgLEf ~ 0.26%. For LSND, the detector is situ- 
ated approximately to the east of the source. This implies 
that the angle between celestial north and the direction 
of propagation of the neutrinos is near 90°, which results 
in the estimate |^| ~ 10^^^. 

In contrast, the KARMEN detector is located roughly 
to the south of the neutrino source, at latitude x — 51°. 
We can therefore approximate Q ~ 180° — x — 129°. 
Taking E = 35 MeV and L = 18 m for KARMEN yields 
a transition probability P,^^^p^ = sin^ 20sin^ 27rL/iyo — 
(sin GgLii^)^ ~ 0.06%. This is more than four times 
smaller than the LSND probability as a consequence of 
the different propagation direction and the smaller dis- 
tance, confirming that direction dependence could help 
reconcile the apparent conflict between KARMEN and 
LSND. 

In the above model, the directional dependence is 

rather limited because the coefficient g introduces only 
@ dependence. This causes minimal variation for any ex- 
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pcrimcnts with both neutrino source and detector fixed 
on the Earth's surface, since the angle 8 is fixed as the 
Earth rotates and is therefore a constant experiment- 
dependent quantity. However, other coefficients can pro- 
duce a strong dependence on $ as well. For instance, 
suppose we choose instead of • The result is an 
oscillation length given by 2-k/Lq = \Esiv? 6cos$^^o|- 
The dependence on $ can substantially change the na- 
ture of an experiment. For purely terrestrial experiments, 
where the source and detector are fixed to the surface 
of the Earth, it fohows that $ = a;©(T - Tq), where 
Cij0 ~ 27r/(23 h 56 min) is the Earth's sidereal frequency 
and To is an appropriately chosen experiment-dependent 
offset. For solar neutrinos, $ varies as the Earth orbits 
the Sun, $ k. Q.q{T - Tq), where = 27r/(l year). 

C. Lorentz-violating seesaw models 

The above models demonstrate some of the strik- 
ing behavior at different energy scales that can arise 
from Lorentz and CPT violation. Mixed energy depen- 
dence among the coefficients for Lorentz violation in /leflF 
can also lead to a Lorentz-violating seesaw mechanism 
that occurs without mass and only in particular energy 
regimes. This can lead to counterintuitive phenomena, 
such as the appearance of a pseudomass in the bicycle 
model of Ref. [29]. In this model, an oscillation length 
emerges at high energies that behaves like a mass-squared 
difference, even though no mass-squared differences exist 
in the theory. 

The bicycle model has nonzero coefficients ^{cLif.^ = 
^{clYJ = 2c and (ai)f^ = (aL)f^ = ajsfi. The ba- 
sic behavior of the oscillation lengths Lah = 2ir/Aab and 
the energy-dependent mixing angle 9 are illustrated in 
Fig. 4. A key feature is that at high energies the line as- 
sociated with the oscillation length L32 resembles that 
from a nonzero mass-squared difference. It turns out 
that the resulting high-energy dynamics reduces to two- 
generation maximal mixing, P^^^^^ ~ sin^(ATO||L/4£^), 
with a Lorentz- and CPT-violating pseudomass Atoq = 
a^cos^ Q/c. 

Unexpected effects of this type can be expected when- 
ever the low- or high-energy limit of hes contains degen- 
eracies. Consider, for example, a 3 x 3 hamiltonian hes 
for which there exists a basis, not necessarily the flavor 
basis, in which we can write 

/ 2hi h2 hsX 
hes= { h*2 , (40) 
\hl J 

where irrelevant diagonal terms are neglected. The 

in teresting eigenvalue difference for this case is A = 
^/ {h\Y + 1/^2^ + |/i3p — hi . Suppose that the mixed en- 
ergy dependence introduced by combinations of masses, 
dimension-one coefficients, and dimensionless coefficients 
enforces hi ^ \/|/i2p + |/i3p at some energy scale. Ex- 




E (GeV) 

FIG. 4: Range of oscillation parameters versus energy in the 
bicycle model with c = 10"^^ and a = 10"^° GeV. (a) Min- 
imum (cos^ 6 = 1) and maximum (cos^ G = 0) of the vari- 
ous oscillation lengths Lab = 2iT/Aab- Note that L32 is un- 
bounded, (b) The allowed range of sin'^ 9 and cos'^ 9 over all 
possible directions, < cos^ 9 < 1, as a function of energy. 



panding the eigenvalue difference then yields [75] A « 

^{\h2\^ + M^)/hi + ---. 

In the bicycle model, /i2 and hs arise from a dimension- 
one coefHcient and are therefore constant with energy, 
but hi arises from a dimensionless coefficient and there- 
fore grows linearly with energy. As a result, at high en- 
ergies the eigenvalue difference is proportional to E~^, 
which resembles the usiial mass case. Using different 
combinations of masses and coefficients for Lorentz vi- 
olation, it is straightforward to construct similar models 
that produce E~^, E~'^, or dependence at high ener- 
gies, or E^, E'^, or E^ dependence at low energies. More 
complicated E^ dependences are possible when the full 
6x6 effective hamiltonian (14) with f <-> i/ mixing is 
considered. 



V. DISCUSSION 

In this paper, we have presented a general framework 
for the study of Lorentz violation in the neutrino sector. 
The key result is Eq. (14), which represents the general 
effective hamiltonian for neutrino propagation in the 
presence of Lorentz and CPT violation. We have ex- 
tracted theoretical implications of this hamiltonian and 
have initiated a study of experimental sensitivities to the 
predicted effects. The various simple models of Sec. IV 
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illustrate some of the key physical features and offer nu- 
merous options for future exploration. 

Our analysis shows that the data from existing and 
near-future neutrino experiments could be used to at- 
tain interesting sensitivities to possible Lorentz- violating 
effects. Moreover, the existing analyses appear insuffi- 
cient to exclude the possibility that some or perhaps even 
all the established neutrino-oscillation signals are due to 
Lorentz violation. 

An interesting open theoretical challenge is to iden- 
tify from the plethora of available choices one or more 
elegant models with features compatible with observed 
data, preferably involving only a small number of degrees 
of freedom. One simple candidate is the bicycle model 
[29] , which has no mass-squared differences and only two 
degrees of freedom rather than the four used in the con- 
ventional massive- neutrino analysis, but which nonethe- 
less reproduces the major observed features of neutrino 
behavior. This and similar models offer one possible path 
to explore, but it is likely that many other qualitatively 
different and interesting cases exist. 

On the experimental front, confirming or disproving 
these ideas would involve analysis of existing and future 
data to seek a 'smoking-gun' signal for Lorentz viola- 
tion. In the remainder of this section, we summarize 
some possible smoking-gun signals and then offer some 
remarks about experimental prospects for detection of 
Lorentz and CPT violation. 



A. Generic predictions 

The numerous options for coefficients for Lorentz and 
CPT violation and the size of unexplored L versus E 
space are impediments to a completely general analysis. 
An alternative strategy to uncover evidence of Lorentz 
violation is to seek model-independent features that rep- 
resent characteristic signals. We list here six classes of 
signal. Confirmed observation of any of them would be 
evidence supporting the existence of Lorentz violation. 

Class I: Spectral anomalies. Each coefficient for Lorentz 
violation introduces energy dependence differing from the 
usual case. Detection of a vacuum oscillation length that 
is constant in E or inversely proportional to E to some 
power would constitute a clear signal of Lorentz viola- 
tion. Note that combinations of masses, dimension-one 
coefficients, dimensionless coefficients, and matter poten- 
tials can produce more complicated energy dependences 
in both oscillation lengths and mixing angles. In general, 
a mixing angle is constant in energy (mly if all relevant 
coefficients for Lorentz violation, masses, and matter ef- 
fects have the same dimension, which requires no more 
than one of these to be present. 

Class II: L-E conflicts. This class of signal refers to any 
null or positive measurement in a region of L-E space 

that conflicts with all scenarios based on mass-squared 
differences. For example, consider a solid line in Fig. 1 



passing through CHOOZ. A measurement of substantial 
oscillation in the D,, sector in any experiment below this 
line would be in direct conflict with a mass-based inter- 
pretation of the CHOOZ results. Signals in this class 
might best be sought by searching for oscillation effects 
in each species of neutrino and antineutrino for regions 
of L-E space in which conventional oscillations are ex- 
cluded. Of the six classes of signal discussed in this sec- 
tion, this is the only one for which there is presently some 
positive evidence, the LSND anomaly. 

Class III: Periodic variations. This class involves signals 
for rotation-invariance violations and contains two sub- 
classes: sidereal variations and annual variations. Con- 
sider first sidereal variations, which have been widely 
adopted as the basis for Lorcntz-violation searches in 
other sectors of the SME. In terrestrial experiments with 
both the detector and the source fixed on the Earth, the 
direction of neutrino propagation relative to the Sun- 
centered frame changes during the sidereal day due to 
the rotation of the Earth. The induced periodic vari- 
ation of observables with time represents a signature 
of Lorentz violation. In the Sun-centered frame, the 
neutrino-propagation angle 6 is constant for a fixed 
source, but the angle $ varies periodically according to 
$ = loq{T — Tq), where Tq is an experiment-dependent 
time at which the detector and source both lie in a plane 
parallel to the XZ plane with the detector at larger val- 
ues of X. The resulting neutrino-oscillation probabili- 
ties exhibit periodic variations at multiples of the side- 
real frequency w®. The second class of periodic signals, 
annual variations, can also arise directly from rotation- 
invariance violation. For solar-neutrino experiments, the 
source is the Sun and the detector changes location with 
time as a consequence of the orbital motion of the Earth 
about the Sun. One can therefore expect variations at 
the Earth orbital frequency il^ and its harmonics. In this 
context, note that the direction p of solar neutrino prop- 
agation in the Sun-centered frame is uniquely given by 
p = (— cos^qT, — cos7?sinn0T, — sinjjsinflQT), where 
r] ~ 23.4° is the angle between the Earth's equatorial and 
orbital planes. We remark in passing that suppressed an- 
nual variations can also arise indirectly as boost- violating 
effects [15, 16, 21, 22] in experiments with terrestrial and 
possibly atmospheric neutrino sources, as a result of the 
noninertial nature of the Earth's motion around the Sun. 

Class IV: Compass asymmetries. This class also results 
from rotation-invariance violations, but the signals are 
independent of time. They can be characterized as the 
observation of unexplained directional asymmetries at 
the location of the detector. For terrestrial and atmo- 
spheric experiments, averaging over time eliminates the 
dependence on the neutrino-propagation angle so the 
result depends only on energy and the angle 0. Rotation- 
symmetry violations can therefore cause a difference in 
observed properties of neutrinos originating from differ- 
ent directions. Note that the east and west directions 
are equivalent under the averaging process, since the $ 
dependence is eliminated, but direct comparison of the 
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north, south, and east directions would be of interest for 
these signals. Note also that the O dependence typically 
introduces vertical up-down effects and could include, 
for example, modifications in the up-down asymmetry 
of atmospheric neutrinos. We remark also that compass 
asymmetries can carry information completely indepen- 
dent of the information in periodic variations. This is 
seen in the example in Sec. IV B 2, which has 8 depen- 
dence but no $ dependence and consequently predicts 
compass asymmetries without sidereal variations. 

Class V: Neutrino-antineutrino mixing. This class of sig- 
nal includes any appearance measurement that can be 
traced to v ^ D oscillation. Any model with nonzero co- 
efficients of type g 01 H exhibits this behavior, including 
the class of simple one-species models discussed in Sec. 
IV Bl. Note that this class of signal involves lepton- 
number violation. 

Class VI: Classic CPT test: Pu.^u^ ^ P^a^Pt- This is 
the traditional test of CPT discussed in Sec. HE, involv- 
ing violation of the result (24a). A related signal would 
be violation of the second result, Eq. (24b), which also 
involves i/ <-> mixing. 



B. Experimental Prospects 

We conclude with some comments about prospects 
for Lorentz- and CPT-violation searches in the major 
types of experiments. Table I provides a summary of the 
present situation. The left-hand part of the table con- 
tains three columns with information about coefficients 
for Lorentz violation. The first column lists combinations 
of coefficient matrices relevant to neutrino propagation, 
extracted from the general hamiltonian (14) and sepa- 
rated according to rotation properties into timelike (T) 
and spacelike (J) components in the Sun-centered frame. 
The second column lists the maximum number of inde- 
pendent degrees of freedom (DOF) associated with each 
combination of coefficient matrices. These numbers can 
be obtained by examining the form of Eq. (14) and us- 
ing the symmetry properties in generation space listed 
in Eq. (4). In certain specific models, some of these de- 
grees of freedom may be unobservable. The third col- 
umn displays the classes of signal that are relevant for 
each coefficient matrix, using the nomenclature of the 
previous subsection. The right-hand part of the table 
contains estimated attainable sensitivities, classified ac- 
cording to each of five types of oscillation experiments. 
Each entry in the table represents the base- 10 logarithm 
of the expected sensitivity to the corresponding coeffi- 
cient for Lorentz violation. The sensitivities shown in the 
table can be obtained by examination of Fig. 1. Given 
an experiment with maximum L coverage of imax and 
maximum E coverage of .Emax! the crude sensitivity a to 
a coefficient for Lorentz violation of dimension 1 — d is 
taken to be cr w — logimax — log -Emax- For simplicity 
in the presentation, it is understood that the sensitivi- 



ties listed for the dimension-one coefficients o,/,, H are 
measured in GeV. The final row of the table contains a 
rough estimate of sensitivities measured in GeV'^^"'*' to 
a generic coefficient {kd)x... for a Lorentz-violating oper- 
ator of nonrenormalizable dimension n = d + 3. Some 
caution is required in interpreting the latter numerical 
estimates because the coefficients {kd)x... are expected 
typically to be suppressed by d-dependent powers of the 
Planck scale. 

The table confirms that Planck-scale sensitivities to 
Lorentz and CPT violation are attainable in all classes 
of experiment, with the most sensitive cases potentially 
rivaling the best tests in other sectors of the SME. Note 
that the estimated sensitivities assume order-one mea- 
surements and therefore may underestimate the true at- 
tainable sensitivity in any specific experiment. Note also 
that a variety of experimental analyses are needed to ex- 
tract complete information on Lorentz and CPT viola- 
tion, with no single class of experiment presently in a 
position to provide measurements of a complete set of 
coefficients. In the remainder of this subsection, we offer 
a few more specific remarks about each type of experi- 
ment. 

Solar-neutrino experim,ents. The abundance and qual- 
ity of the current solar-neutrino data make these experi- 
ments a promising avenue for Lorentz- violation searches. 
The relatively large range of solar-neutrino energies sug- 
gests interesting information about spectral anomalies 
might be obtained, but complications introduced by mat- 
ter effects are likely to make this practical only in rela- 
tively simple cases such as the FC model (28). Of the 
other classes of signals, periodic variations and neutrino- 
antineutrino mixing may be the most relevant to solar 
neutrinos. The periodic variations in observables would 
occur at multiples of Jl^, appearing despite compensa- 
tion for the flux variation due to the eccentricity of the 
Earth's orbit. Direct detection of any antineutrinos orig- 
inating from the Sun would be evidence of •<-> mixing 
and hence of possible Lorentz violation. 

Atmospheric-neutrino experiments. Like solar neutrinos, 
atmospheric neutrinos cover a relatively large region of 
L-E space, but complications from matter effects hin- 
der a general spectral-anomaly search. However, Fig. 2 
shows that searches for atmospheric oscillations at the 
highest energies and largest distances could reveal oscil- 
lations absent in the usual solution, thereby providing 
evidence for L E conflicts. Atmospheric neutrinos orig- 
inate from all directions, so they are an ideal system for 
directional-dependence searches. Not only are they sen- 
sitive to sidereal variations, but also the directional capa- 
bilities of detectors like SK make atmospheric neutrinos 
perhaps the most promising place to search for compass 
asymmetries. Moreover, since atmospheric data involve 
both neutrinos and antineutrinos of two species in compa- 
rable numbers, it may be possible to address both v ■(-^ i> 
mixing and the classic CPT tests (24a) and (24b). 

Reactor experiments. Nuclear reactors are good sources 
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of Vei and they are therefore well suited to searches for 
y 9 mixing. Since both the sources and the detectors 
in all these cases are fixed, the experiments are also sensi- 
tive to sidereal variations, and some may have additional 
sensitivity to compass asymmetries. For example, the 
reactor experiment KamLAND detected neutrinos from 
multiple reactors and different locations. Experiments 
with multiple sources like this can analyze their data for 
compass asymmetries that depend on the direction to the 
various neutrino sources. 

Short-has eline accelerator experiments. LSND already 
seems to suggest a positive L-E conflict, which will be 
tested by the forthcoming results of the MiniBooNE ex- 
periment. Many of these short-baseline accelerator ex- 
periments are especially interesting for signals based on 
L-E conflicts because they operate in a region of L-E 
space where the conventional mass scenario predicts no 
oscillations. Sidereal variations can readily be sought by 
experiments such as CHORUS, KARMEN, MiniBooNE, 
NOMAD, and NuTeV, since each has a fixed source and 
detector. Note that the existing data from these exper- 
iments could in principle contain a positive signal for 
sidereal variations because the published null results arc 
based on an average over time. The well-defined flavor 
content of the sources for these experiments may also of- 
fer sensitivity to v v signals and to the classic CPT 
test. Sonic of these experiments, such as MiniBooNE 
and NuTeV, may be particularly sensitive to Lorentz vi- 
olation because they can switch from a predominately 
source to a predominately i?^ source. 

Long-baseline accelerator experiments. Several future 
long-baseline accelerator-based experiments, such as 
ICARUS, MINOS, and OPERA, are planned to probe 
the GcV region oi L E space at distances of hundreds of 
kilometers, and some results in this regime have already 
been reported by K2K. These experiments can search 
for oscillations in obtained from meson decays, and 
they are designed to test the atmospheric-oscillation hy- 



pothesis. Nonetheless, L-E conflicts are still possible: 
a measurement of — » z/g, for example, would repre- 
sent an L-E conflict because this oscillation is absent at 
these energies and distances in the conventional scenario 
with masses. The data obtained can be also analysed 
for sidereal variations, since in each case the source and 
detector arc fixed. Moreover, except for OPERA and 
ICARUS, which are both part of the CERN Neutrinos to 
Gran Sasso (CNGS) project, the beamline for each ex- 
periment points in a different direction. This means each 
is expected to respond differently to rotation-invariance 
violations. These experiments may also be able to ad- 
dress V y mixing and the classic CPT signal, since the 
flavor content of the beams is well known. 

Other experiments. Experiments designed to search for 
neutrino properties other than oscillations can also ad- 
dress Lorentz violation. To some extent, most experi- 
ments are sensitive to sidereal variations and compass 
asymmetries. The other signals discussed in Sec. VA 
are more unique to neutrino oscillations, but analogous 
signatures arc likely to arise in most cases. 

One possible test of Lorentz invariance involves a di- 
rect comparison of velocities of supernova neutrinos and 
photons, such as those from SN1987A [34, 76], which 
could be performed either by some of the experiments 
listed above or by neutrino telescopes. A similar method 
has been applied in the photon sector, where the veloci- 
ties of different polarizations arc compared [22] . Another 
method that could be adapted to the neutrino case is a 
simple pulse-dispersion analysis. The energy dependence 
and the independent propagation of each he« eigenstate 
imply that different components of the neutrino pulse 
propagate at different velocities, causing the pulse to 
spread. For SN1987A, all the observed neutrinos arrived 
in a time interval of about 5T ~ 10 s and had energies 
E :^ 10 — 20 MeV. Since these neutrinos took roughly 
To ~ 5 X 10^^ s to reach the Earth, we can crudely esti- 
mate that the maximum difference in velocity across the 
AE ~ 10 MeV energy spread of the /leff eigenstates is 
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Svjc ~ 5T/Tq ~ 2 X 10~^^. Wc can then make a simple 
dimensional estimate of the sensitivity of this method to 
various terms in /leff • This suggests a sensitivity of about 
^200 eV to mass terms, 2 x 10~^^ GeV to dimension-one 
coefficients, and 2 x 10~^^ to dimensionless coefficients. 
The mass estimate agrees with the result of a detailed 
analysis along these lines [77]. 

Lorentz violation may also be relevant to direct mass 
searches such as the proposed KATRIN experiment [32], 
designed to measure directly the mass to better than 
1 eV. Within the currently accepted solution to the os- 
cillation data, a mass matrix with eV-scale masses but 
mass-squared differences of 10~^ eV^ and 10~^ eV^ would 
be nearly degenerate. This seems unlikely in light of the 
charged-lepton mass hierarchies. However, siippose that 
the mass matrix is nearly diagonal and that neutrino os- 
cillations are primarily or entirely due to Lorentz vio- 
lation instead. Then, individual masses of eV order or 
greater may be present with little or no effect on the ex- 
isting neutrino-oscillation data, but they would produce 
a signal in the KATRIN experiment. 

Another area of widespread interest is the search for 
neutrinoless double-beta decay [33]. This decay mode is 
an indicator of lepton-number violation, which can result 
from Majorana-type couplings introduced by Majorana 
masses or by gauge- violating coefficients for Lorentz vio- 
lation. Many of the null results of searches for neutrino- 
less double-beta decay could therefore be reanalysed to 
yield constraints on certain types of Lorentz violation. 
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APPENDIX A: EFFECTIVE HAMILTONIAN 



This appendix presents some details for the deriva- 
tion of the effective hamiltonian (14). We first perform 
a spinor decomposition of the hamiltonian in the mass- 
diagonal Majorana basis. The result is then block di- 
agonalized in the light-neutrino sector and transformed 
into the original weak- interaction basis. We remind the 
reader that generation indices in the mass-diagonal basis 
are A' = 1, . . . , 6 for = 3 neutrino species, while the 
restriction to light neutrinos in this basis is represented 
by indices a' = 1,2,3. Also, in the flavor basis, upper- 
case indices take the values A = e,iJ,, r, e*^, /x*-^, r*^, while 
lower-case ones span a = e,i2, r. 



1. Spinor decomposition 

In this section, we project the hamiltonian onto the 

massless spinor basis used in Eq. (10). This corresponds 
to choosing a convenient p^dependent 7-matrix basis that 
allows us to write the equations of motion in terms of the 
b and d amplitudes. 

Working in the mass-diagonal basis, the hamiltonian 
is given by 

Ua'B'{p) = U<iA'B'{p) + 5Ua'B'{p), 

'HoA'B'if) = 7°(7-P + 'Tl(A'))^A'S' , 

6Ha'b'{p) = -^{7°Sr°Hoip) +Ho{p)7°Sr°)A'B' 



+7°(5f • p + 6M)a'b' 



(Al) 



It turns out to be useful to decompose r^,g, and Ma'B' 
in terms of 7 matrices, as in Eq. (2). Therefore, we write 

r:;,^, = ^°UA'Ai°r'Xs{UB'B)^ 

= YSa'B' + c'^A'B'l^i + d'X^,B,'j5lf, 
+^A'B' +«/^'B'75 + 29A'B''^^li ' 

Ma'b' = 7°Ua'aJ°Mab{Ub'b)^ 

= niA'B' + im5A'B'^5 

+«A'S'7^ + &^i's'757^ + hK'B'^^f^-^ ■ (A2) 

We begin the spinor decomposition of the hamiltonian 
(Al) by considering the Lorentz-covariant terms. The 
properties of the massless spinor basis imply that the 
only nonzero projections of the kinetic term are 

UL,R{'P)b^1 ■ P ^A'B')ul,r{'P) 

= -'"r,l i-P) (7°7 • P Sa'B' )vr,l {-p) 

= \p\5a'B' , (A3) 

while the surviving projections of the mass term are 
"L,fl(P) {'y°'m{A')SA'B' )vr,l {-P) 

= UL,R{p)VR,L{-p)'m(^A')^A'B' (A4) 

and conjugates. The quantities ul,r{p}''^b..l{^P) are 
phases that can be chosen arbitrarily by changing the rel- 
ative phase between ul,r{p) and Vr^l^—P) = f^i—p)- 
For the spinor decomposition of the Lorentz-violating 
terms in the hamiltonian (Al), we define the 2x2 ma- 
trices 

AA'B'{p)=A^B'A'iP) 

_ f ' 



u^r(p) 



6nA'B>{p){uL{p),UR{p)), (A5) 



Aa'B'{p) = -A.I,a'{-P) 
= (^"^^^y-HA'B'ipjivRi-pj^vU-p)). (A6) 

It can be shown that the mass-basis analogues of the 
relations (3) are r^,^, = -C{Vg,^,fC-^ and Ma'B' = 
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C{Mbi A')^ C^^ . Note that this corresponds to C ^ /, 
which reflects the Majorana nature of neutrinos in this 
basis. These identities may then be used to show that 
C^l°nA'B'{p)l''C = -[HA'B'i-p)]*. Finally, with the 
aid of the relation va^Lip) = CuJ^ r{p), it follows that the 
remaining terms in the spinor decomposition are given in 
terms of A and A by 



-A 




5nA'B'{p){uL{p),UH{p)), (A7) 



SnA'B'{p){vR{-p),VL{-p)). (A8) 



This implies that the 2x2 matrices Aa' b' , A^/ b' deter- 
mine the Lorentz- violating effects. 

Combining the above results, we obtain the spinor- 
decomposed hamiltonian appearing in Eq. (12): 



where rj is the 2x2 diagonal matrix of phases T}{p) = 
-T]{-p) = diag[uL{p)vR{-p),UR{p)vL{-p)]. 



We seek an explicit expression for Aa'b'- The next 

subsection shows that the effects of Aa'B' are sub- 
leading order, so we concentrate here on the projec- 
tions in Aa'B', which involve the spinors Ul and ur. 
It is useful first to find expressions for the quantities 
^'a{l)75)7^)757'')Cr'"'}^t/3) where a, /3 = L,R. We ob- 
tain the following nonzero results: 



P'"5al3l\p\, 

{uRa^^'^ULf 

is/2{j^{e+r-p''{e+Y)/\p\, 



Ul(J^"Ur 



(AlO) 




Ia'{p) 

\P\ 'rn(^A')V{p) 
■mi^A')V*{-P) -\p\ 



A'B 



A 



A'B 



iP) 



A'B 



'i-P) -K'B'i-P) 



(A9) 



where Sl = 1, Sr = —1, p^ = {\'p\;p), and (e+)'' satisfies 
the relations (15). With these results and Eqs. (Al) and 
(A2), we can extract the projections of 6H onto ul and 
Ur- 



Aa'B' = 




(a + b^pf, - (c + dy-'p^p.U'B' -iV2pf,ie+)4g^-'^pa - H'^-'U'B' 
i^/2p^{e+)l\g^"''p, - H^'^Ub' [{a - hYp^ - (c - dyp^^P.U'B' 



(All) 



r 



In this expression, we neglect off-diagonal terms entering 
as mass multiplied by coefficients for Lorentz violation, 
since in most situations these terms are suppressed rela- 
tive to those above. 



2. Block diagonalization 

The above spinor decomposition of the hamiltonian is 
independent of the specific neutrino mass spectrum. To 
make further progress, we adopt the scenario described 
at the beginning of Section II B and restrict attention to 
ultrarelativistic dynamics in the subspace of light neutri- 
nos, spanned by the a' indices. The hamiltonian is then 
dominated by the diagonal kinetic term in Eq. (A9). The 
upper and lower diagonal blocks of this term have oppo- 
site sign, so they differ by an amount large compared 
to both mass and Lorentz- violating terms. This in turn 
implies that standard perturbation techniques to remove 
the off-diagonal blocks can be applied. As a result, terms 
in the off-diagonal blocks of Eq. (A9) appear at second 



order in the block-diagonalized form. One consequence is 
that the leading-order mass contribution appears at sec- 
ond order, whereas certain forms of Lorentz violation ap- 
pear already at first order. This feature can ultimately be 
traced to the 7-matrix structure of the Lorentz-covariant 
portion of the theory. 

Provided the conditions TO(a'), |Aa'6'||Ao/h/| <C \p\ are 
satisfied, the block diagonalization of Eq. (A9) can pro- 
ceed through the perturbative construction of an appro- 
priate unitary matrix U. First, write lA in the form 
U — I + e^"'^-* -|- e^^-* -|- . . ., where e^") is of nth order in the 
dimensionless small quantities m(^a')/\p\, Aa'b'/\p\, and 
^a'b' /\p\- The block-diagonal hamiltonian resulting from 
this transformation can be expanded in a similar fashion: 



a'b' 



Ud'a'Ha'h'lAyy 



if (0) + H"-'' 

ah ah 



r(i) 



ah 



(A12) 



where each H^""- is nth order in small quantities. The 

a'h' 

zeroth-order term H^'^l is the usual kinetic term, which 

a'b' 
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is already block diagonal. The first-order term H^]^, can matrix Va'a, we find 



be made block diagonal by an appropriate choice of e^^^. 
A suitable leading-order transformation is 



where 



(1) ^ Sa>a> I Ca'b' (p) 



ea'b'ip) = m(^a')Sa'b''ri{p) + K'b'ip)- 



(A13) 



(A14) 



Using e^^^ and -ff^f^,, which depends on both e'^^ and e'^', 

we can find e^^^ and then continue iteratively to arbitrary 
order. 

Under the transformation U, the hamiltonian re- 
stricted to light neutrinos may be written 



a'b' 







-hZ,y{-p) 



(A15) 



Calculating U to second order in small quantities yields 
the second-order hamiltonian 



(A16) 

This expression neglects terms that are second order in 
coefficients for Lorcntz violation and terms that enter as 
the product of TO(a')/|p| with A. The latter terms con- 
stitute subleading-order corrections under the reasonable 
assiimption that A and A arc comparable in size. 

While formally the two bases related hyU are different, 
in practice this difference is of little consequence. Our 
main goal is to determine oscillation probabilities. The 
effects of lA appear in the mixing matrix and therefore 
modify the amplitudes of oscillations. However, since U 
is close to the identity, the basis change produces only 
tiny and unobservable changes in oscillation amplitudes. 
It therefore suffices in practice to assume U = I for pur- 
poses of the basis transformation, corresponding to ignor- 
ing the difference between the a' and a' indices. Similar 
arguments apply to the field redefinition relating v and 
X- This also underlies the validity of assuming unitarity 
mixing matrices in the conventional case with neutrino 
mass, even though the siibmatrix Va'a is only approxi- 
mately unitary. In contrast, the diagonalization of h in 
Eq. (A16) can introduce arbitrary amounts of mixing. 

The above description in the mass-diagonal basis com- 
pletely determines the neutrino dynamics, but in practi- 
cal situations a description in the weak-interaction basis 
is more useful. This requires the transformation of ha'b' 
to the original flavor basis. 

The first step in implementing the desired transforma- 
tion is to determine the relation between the coefficients 
in Eq. (2) and those in Eq. (A2). In terms of the unitary 



"-A'B' 

ifl'B' 
1 „Vi' 

rUA'B' 

imsA'B' 

C-A'B' 

"A'B' 
2^ A'B' 



ReVA'AV^,B{c + dr^j,, 
ilm Va'aV|,b(c + d)%, 

ilmVA'AVB'B[{e + ifYC]AB, 
ReVA'AVB'B[{e + l.fYC\AB, 
^eVA'AVB'B\{g^^''C)AB 

-lmVA'AVB'Bke^'""'(9p'T''C)AB, 
BeVA'AVB'B[{m + im5)C]AB = m(^A')SA'B', 
ilmVA' aVb' B[{m + i'm^)C]AB = 0, 
ilm VA'AVB'Bia + &)as> 
ReVA'AV^'Bia + bY^s, 

llmVA'AVB'B\{H^'''C)AB 

+iReVA'AVB'By'"'''^iHpX)AB. (A17) 



Note that all the coefficients in the mass-diagonal ba- 
sis are either pure real or pure imaginary, reflecting the 
Majorana nature of neutrinos in this basis. Using this 
equation, we obtain 

[{a+brp^-{c+drypP,]a'b' 

= [{a + bYp^ - (C + dY''p^.pAabVa'aV:,^ , 

[{a-bYPf.-{c-dY''p^pAa'b' 

= [-{a + bYp^ - (c + dY-p^^PujX^Vb'b , 
-iV2p^{e+Yrp,-H^^''U' 

= -is/2p^{e+Y[{g>""'p„ - H>"')CUVa'aVb'b , 
iV2p^{e+YA9'""'p. - H^wb' 

= i^np^ie^YAir^P. + ancMbKa^b , (A18) 

using the assumption that the submatrix Va'a is unitary. 

Within a standard seesaw mechanism, the right- 
handed Majorana-mass matrix R appearing in Eq. (8) 
is large, \R\ » |i|,|£>|. Calculating the matrix Vab at 
leading order in small mass ratios and |Z)|/|i?| 

produces the identity 

m(^a')Sa'b' = Va'aVb'b{mi)ab, (A19) 

where mi = L — DR~^D^ , and hence the relation 

m'l^')5a'b' = Va'aVC'b{'mim\)ab = V*,^Vb'b{mim\Yab- 
Combining results yields the desired form. 



(A20) 



1^ + 2^'^^))'^"''''+^''''''^^^ 



Va'a 

v*„ 



{hes)ab 



V* 

^b'b 




Vb'b 



(A21) 



where h^s is given in Eq. (14). 
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APPENDIX B: MINIMAL SME TERMS 

Restricting attention to the coefficients (cl)^^, {aL)ab^ 
which arc contained in the minimal gauge-invariant SME, 
effectively decouples neutrinos and antineutrinos and 
produces vanishing transition probabilities (21c) and 
(21d). This appendix describes a useful parametrization 
of these coefficients. 

Each coefficient matrix for Lorentz violation can be 
parametrized with three eigenvalues and a constant uni- 
tary matrix. We define 

\ 

{cLr = (Un^ {cl)^^^ U^'^ (Bl) 

V icr^^J 

for each coefficient matrix {clY" , and 

\ 

[aLT = {Un^ (aL)f2) f/" (B2) 

V {aL)'f,J 



for each coefficient matrix {a^)^. The unitary diagonal- 
izing matrices U^^ , are chosen so that if there is only 
a single nonvanishing coefficient matrix then UeS in Eq. 
(19) takes the block-diagonal form 

The reader is warned that the above decomposition is 
frame dependent, so neither the eigenvalues nor the mix- 
ing matrices behave as tensors under observer Lorentz 
transformations. We therefore advocate restricting this 
type of decomposition to the standard Sun-centered ce- 
lestial equatorial frame. 

Adopting a CKM-like decomposition of the U matrices, 

we denote mixing angles and phases associated with each 
{clY" by e^^:^, 0^'3^ 0^3^ and /^r, 13^"- Similarly, for 
each {aLY we write 61^2, 6^^^, 6^^, and S"", /3f , /3^. The 
U matrices may then be written explicitly in the form 



■-12 ^-13 
*12 '^13 
*13 ^ 



^12 '--13 
*12'^13 



■^12 '--23 
''12 ''23 



'12^23 



'12 ''23 ''13 ' 
^ly 

^12 -^23 '^13 ^ 
^23 ''13 



„M p-iiS" 



12 '-23 



'12^23 ''13 
M M M 
■Sl2'S23'5l3 

'23 ''13 



'12 ■'23 
''12 *23 



^12 ''23 '^13 ^ 

''12 ''23 •'13 ^ 
''23 ''13 



■'12''23 
''12''23 



' '^12'^23'^13^ 
''12 ''23 ''13'^ 
''23 ''13 



1 

e^^r 





1 






























(B4) 



(B5) 



where s^^' = sin6'^^^ c^^" = cos6';;^^ s^^ = sin6l^;„ and 
c^6 = cosC,. 

In the conventional massive-neutrino analysis, the /? 
matrix of phases can be absorbed into the amplitudes 
ba{t;p) and da{t;p), so these phases are normally unob- 
servable and can be neglected. However, in the present 
context, only one set of /3 phases may be absorbed in this 
fashion. The presence of multiple coefficient matrices for 
Lorentz violation implies that they cannot typically be 
neglected. 

Neutrino oscillations arc insensitive to terms in the ef- 
fective hamiltonian that are proportional to the identity. 
Consequently, only two eigenvalue differences for each co- 
efficient matrix for Lorentz violation contribute to oscil- 



lation effects. Also, each coefficient matrix is associated 
with three mixing angles and three phases. It follows 
that the maximum number of gauge-invariant degrees of 
freedom that enter into neutrino oscillations in the min- 
imal SME alone is 16x8 for cl and 4x8 for a^, for a 
total of 160. However, some of these are unobservable. 
The 8 trace components ?7pjy(ci)^'^ are Lorentz invariant, 
and both these and the 6 x 8-component antisymmetric 
piece of {clY" are absent in the leading-order hamilto- 
nian (14). This leaves 104 leading-order degrees of free- 
dom in fli and c^, in agreement with the numbers listed 
in Table I. For the minimal SME, one set of /3 phases 
is also unobservable, which reduces the total number of 
degrees of freedom to 102. 
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